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Abstract 

We derive semiclassical asymptotics for the orthogonal polynomials Pn{z) 
on the line with respect to the exponential weight exp(— A^y(z)), where V{z) 
is a double-well quartic polynomial, in the limit when n,N ^ oo. We assume 
that e < {n/N) < Acr — e for some e > 0, where Acr is the critical value which 
separates orthogonal polynomials with two cuts from the ones with one cut. 
Simultaneously we derive semiclassical asymptotics for the recursive coefficients 
of the orthogonal polynomials, and we show that these coefficients form a 
cycle of period two which drifts slowly with the change of the ratio n/N. The 
proof of the semiclassical asymptotics is based on the methods of the theory 
of integrable systems and on the analysis of the appropriate matrix Riemann- 
Hilbert problem. As an application of the semiclassical asymptotics of the 
orthogonal polynomials, we prove the universality of the local distribution of 
eigenvalues in the matrix model with the double-well quartic interaction in the 
presence of two cuts. 
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1. Introduction and formulation of the main theorem 

About thirty years ago Freeman Dyson found an exact solution for the 
scaling limit of correlations between eigenvalues in the Gaussian unitary en- 
semble of random matrices. He conjectured that this scaling limit should 
appear in a much broader class of non-Gaussian unitary ensembles of random 
matrices. This constitutes the famous conjecture of universality in the theory 
of random matrices. Dyson found also a remarkable formula which expresses 
the eigenvalue correlations for finite N (before the scaling limit) in terms of 
orthogonal polynomials on the line with respect to an exponential weight (see 
[Dys]; see also [Meh], [BIPZ], and others). This reduces the universality con- 
jecture to semiclassical asymptotics of orthogonal polynomials. 

The theory of orthogonal polynomials is a classical branch of mathemat- 
ical analysis and it finds applications in different areas of pure and applied 
mathematics. Among the most exciting recent applications of orthogonal poly- 
nomials related to the random matrices are those to quantum gravity, string 
theory, and integrable PDEs (see, e.g., [ASM], [Dem], [Wit] and references 
therein). In the present paper we show how these recent developments can be 
beneficial for the theory of orthogonal polynomials itself and, in particular, for 
the aspects of the theory related to the matrix models. 

The Dyson formula gives the eigenvalue correlations for finite A'^. To get 
the universal correlation functions in the limit when N ^ oo one needs to 
obtain semiclassical asymptotics for orthogonal polynomials. For the Gauss- 
ian matrix model the problem reduces to the semiclassical asymptotics for the 
Hermite polynomials which were obtained in the classical work of Plancherel 
and Rotach [PR]. The Hermite polynomials are the eigenfunctions of the quan- 
tum linear oscillator (up to a Gaussian function) and their asymptotics follow 
from the semiclassical formulae for the Schrodinger operator with quadratic 
potential. Using the Plancherel-Rotach asymptotics in the inner region, be- 
tween the turning points, Dyson evaluated the scaling limit of correlations in 
the Gaussian matrix model inside the support of the limiting spectral measure. 
The scaling limit is expressed in terms of the sine-kernel. Dyson's conjecture 
of universality states that the same scaling limit in the bulk of the spectrum 
should appear in the general case of non-Gaussian unitary ensembles of ran- 
dom matrices. The conjecture of universality was then extended to the scaling 
limit of correlations at the edges of the support of the limiting spectral mea- 
sure. For the Gaussian matrix model the scaling limit of correlations at the 
edge is expressed in terms of the Airy kernel (see [BB], [For], [TW2]), and 
this follows from the Plancherel-Rotach asymptotics near the turning points. 
The conjecture of universality states that the same scaling limit is to appear 
at all regular edges of the support of the limiting spectral measure for a gen- 
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eral unitary matrix model. Here regular means that the spectral density has 
a square-root singularity at the edge. The conjecture of universality has, in 
fact, an even more general formulation. One should expect universal scaling 
limits to exist at regular critical points, tricritical points, etc., each being re- 
lated to its own class of universality. The critical points do not occur in the 
Gaussian matrix model but they do in non-Gaussian matrix models, and the 
critical points correspond usually to bifurcations of the support of the limiting 
spectral measure. In this paper we prove the universality conjecture inside 
the spectrum and at the edges for the matrix model with quartic interaction. 
We will discuss the universality at the critical points in another publication. 
A very different approach using estimates for resolvents, to the universality 
conjecture for a large class of potentials has been given in the paper [PS] by 
Pastur and Shcherbina. 

We derive the scaling limit of correlations from the semiclassical asymp- 
totics of orthogonal polynomials. For the general matrix model the problem of 
semiclassical asymptotics is much more difficult than for the Gaussian model. 
The corresponding Schrodinger operator is very complicated and, what is even 
more important, it depends on unknown recursive coefficients for orthogonal 
polynomials. So one has to solve the two problems simultaneously: to find 
semiclassical asymptotics for the recursive coefficients and for the orthogonal 
polynomials. The recursive coefficients for orthogonal polynomials satisfy a 
nontrivial nonlinear equation which was discovered independently in mathe- 
matics and in physics. Mathematicians call it the Freud equation (and its 
generalizations), in the name of Geza Freud who first derived it for orthogonal 
polynomials with some exponential weights on the line. For physicists this 
equation is known as a discrete string equation, and it was derived and studied 
in the papers by Bessis, Brezin, Itzykson, Parisi, and Zuber, and others (see, 
e.g., [BPIZ], [BIZ], [IZ]), devoted to the problem of enumeration of Feynman 
graphs in string theory. The Freud equation does not admit a direct solution, 
and this is the main reason why despite many interesting and deep achieve- 
ments in the theory of orthogonal polynomials (see, e.g., [LS], [Mag4], [Nev3] 
and references therein), the semiclassical asymptotics for general exponential 
weights remains one of the central unsolved problems in the theory. Our new 
technique is based on the important observation made in [FIKl] that the Freud 
equation is integrablc in the sense that it admits a Lax pair. This motivates 
use of the powerful tools of the theory of integrable systems, the method de- 
veloped in the present paper. The central role in our scheme is played by the 
Riemann-Hilbert setting of orthogonal polynomials suggested in [FIK4]. 

To formulate the problem let us consider a polynomial V{z) with real 
coefficients, 

V{z) = ao + aiz H h 02^2:^^, a2p > 0, 
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and orthogonal polynomials on the real line, 

P„(z) = z" + ..., n = 0,1,2,..., 
with respect to the exponential weight e~^'^^^^dz, i.e.. 



/•oo 



(1.1) / P„(z)P„(z)e-^^(^)ci2 = M, 



hn > 0. 



J —oo 



We normalize the orthogonal polynomials by taking the leading coefficient 
equal to 1. The number N in the weight is a large parameter and we are inter- 
ested in the asymptotics of the polynomials Pn{z) in the limit when n, AT — > oo 
in such a way that 

n 

for some fixed e > 0. In this approach the large parameter N plays the same 
role as the inverse Planck constant plays in the semiclassical asymptotics of 
quantum mechanics. 

The orthogonal polynomials Pn{z) satisfy the basic recursive equation 



where Qn and Rn are some coefficients which depend on N. The semiclassical 
asymptotics of the coefficients Q„ and Rn are tied to one of the polynomials 
Pn{z) and we solve these two problems together. 

Although our approach is quite general and can be applied to a general 
polynomial V{z), in this paper we will work out all the details for the quartic 
polynomial 



and we will consider the most interesting case when t < {a double-well 
potential). Since V{z) is even, (1.2) is simplified to 



(1.2) 



zPn{z) = Pn+l{z) + QnPn{z) + RnPn-l{z) 



(1.3) 



g>o, 



(1.4) 



zPn{z) = Pn+l{z) + RnPn-l{z), 



where 



(1.5) 




In addition, integration by parts gives 

(1.6) P^(2;) = NRn[t + g{Rn-i + P„ + P„+l)]Pn-l(z) 

+ g{NRn_2Rn-lRn)Pn-3{z), = ^ • 

Since Pn{z) = nz'^~^ + . . . , this implies the Freud equation 

(1.7) n = NRn[t + g{Rn-i + Rn + Rn+i)l n>l, 
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(cf. [Fre]). The difference equation (1.7) is supplemented by the foUowing 
initial conditions: 



-Ro — 0, i?i 



Prom (1.5) and (1.7) it follows that 



(1.8) 

Let 
(1.9) 

Then 
(1.10) 



0<Rn< 



25 



A- - 



Pn{z)e 



-NV{z)/2 



Our main goal is to prove semiclassical asymptotics for the functions il)n{z) 
and for the coefficients Rn in the limit when A^, n ^ oo in such a way that 
there exists £ > such that the ratio X = n/N satisfies the inequalities 



(1.11) 

where 
(1.12) 



£ < A < Acr — £, 



A = — 



n 

iV' 



In what follows the potential function 

Uiz)-~^^^'''^'^' 



is important. We introduce the turning points zi and Z2 as zeros of U{z) 



(1.13) 



^1,2 



The condition (1.11) implies that zi and Z2 are real, and Z2> z\> C^fe. We 
prove the following main theorem. 

Theorem 1.1. Assvme that N,n oo in such a way that (1-11) holds. 
Then there exists C = C{e) > such that 



(1.14) 



Rn — 



(-l)Vt2-4Ag 
2ff 



< CN-\ A 



n 
iV 
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In addition, for every S > 0, in the interval zi + 5 < z < Z2 — S, 
(1.15) 



where 

and 
(1.16) 



^/sm(f) 



COS 



n 



+ 5) fsm2(p 



2 V 2 
(f) = arccos q, % = arccos r, 



(-1)"X vr 
4 4 



+ o(iv-i) , 



gz'+t 



2y/Yg-tq 



2./Yg ' 2./Ygq-t 
If z > Z2 + S or < z < zi — 6, then 
(1.17) 



V'n(2) = (-1) 



"v/ sinh (f) 



exp < - 



(n+i) 
2 

(-i)"x 



A' 



sinh(2(;!!>) 



n + ; 



+ 



N{l + \z\) 



where 



a = < 



Z > Z2 + S, 

'W 

L 12. 



A; if Q < z < zi — 5 and n = 2k or 2A; + 1, 



cosh \q\, X = cosh |r| 



where q,r are given by (1.16). 

If Zj — S < z < Zj + S, j = 1, 2, i/ien 
(1.18) 



Mz) = ^4^, f(l + n(^)) Ai (iV2/3u;(2)) + r2(z) Ai' (n^/'w{z)) 
V\w'{z)\ L V / \ / 

where ri{z) = O (^^^) , ?'2(-z) = O (A^^^/'^), Ai (z) is t/ie Airy function, and 



W[Z] IS 



(1.19) 



an analytic function on [zj — 6, zj + 6] such that for {z — Zj^^){—iy > 0, 

2/3 



W{Zj 



\/Un(v) dv 

where z^ = Zj + 0{N~^) is the closest to zj zero of the polynomial 
(1.20) Un{z) = U{z) + iV-i (^-^ + ^ + 



= z 



'{gz'^ + t)' 
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and 



-/-(-l)'V/^-4Ag 
25 



The constant factor Cn in (1.15) and (1.17) satisfies the asymptotic equation 
and Dn in (1.18) is 

(1.22) D„ = ivV6^(-l)-o, (7o = (2-j) 

Finally, the asymptotics of hn in (1.1), (1-3) is 



n 
L2 



(1.23) 



hn = 27r-\/^ exp 



Nt'^ NX 



(1 + In I) +0(iV-i) 



The asymptotic formulae (1.15), (1.17), and (1.18) are an extension of the 

classical Planchcrcl-Rotach asymptotics for the Hermite polynomials (see [PR] 
and [Sze]), to the orthogonal polynomials with respect to the weight e"^^*^^-* 
where V{z) is the quartic polynomial (1.3). These formulae are extended into 
the complex plane in 2; as well (see Theorems 7.5, 7.7 in Section 7 below). The 
formula (1.15) can be rewritten as the semiclassical formula 



(1.24) M^) 



|C/(z)|V4 



cos 



^ I \UNiv)\^^^dv + ^ I +0{N 



where Un{z) is as defined in (1.20). Asymptotics (1.14) of the coefficients Rn 
are Preud-type asymptotics. For the homogeneous function V{z) = and 
some of its generalizations, the asymptotics of Rn are obtained in the papers 
of Freud [Pre], Nevai [Nevl], Magnus [Magi, 2], Lew and Quarles [LQ], Mate, 
Nevai, and Zaslavsky [MNZ], Bauldry, Mate, and Nevai [BMN]. Semiclassical 
asymptotics of the functions ij^niz) are proved then for V{z) = z'^ by Nevai 
[Nevl] and for V{z) = by Sheen [She]. Por general homogeneous V{z), 
somewhat weaker asymptotics arc obtained in the works of Lubinsky and Saff 
[LS], Lubinsky, Mhaskar, and Saff [LMS], Lubinsky [Lub], Levin and Lubinsky 
[LL], Rahmanov [Rah], and others. Application of these asymptotics to ran- 
dom matrices is discussed in the work of Pastur [Pas] . The distribution of zeros 
and related problems for orthogonal polynomials corresponding to general ho- 
mogeneous V(z) are studied in the recent work [DKM] by Deift, Kriecherbauer, 
and McLaughlin. Many results and references on the asymptotics of orthogo- 
nal polynomials are given in the comprehensive review article [Nev3] of Nevai. 
The problem of finding asymptotics of i?„ for a quartic nonconvex polynomial 
is discussed in [Nev2, 3], and is known as "Nevai's problem." 
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Equation (1.14) shows that if < A < Acr then 

-t - ^yt^ - AXg 



(1.25) hm i?2m = L{X) : 

hm -R2m+i = -R(A) 

iV^oo; (2m+l)/Ar^A 



25 



-t + _ 4_xg 
25 



Both L{\) and i?(A) satisfy the quadratic equation 

gu^ + tu + \ = 0, 

so that, when n grows, i?„ jumps back and forth from one sheet of the parabola 
to another (see Fig.l). In other words, the sequence {Rn-, n = 0,1,2,...} 
forms a period two cycle which is slowly drifting with the change of the ratio 
n/N. At A = Acr the two sheets of the parabola merge, i.e., L(Acr) = -R(Acr). 
For A > Acr, 



(1.26) 



where u = Q{X) satisfies the quadratic equation 

3gu^ + tu- X = 




Figure 1. The qualitative behavior of the recurrence coefficients 
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(which foUows from the Freud equation (1.7) if we put u = Rn-i = Rn = 
Rn+i)- We consider semiclassical asymptotics for A > Acr and in the vicinity of 
Acr (double scaUng hmit) in a separate work. The difference in the asymptotics 
between the cases A < Acr and A > A^ is that for A < Acr the function ipniz) 
is concentrated on two intervals, or two cuts, [—Z2,—zi] and [zi,Z2], and it is 
exponentially small outside of these intervals, while for A > Acr, zi becomes 
pure imaginary, and ipniz) is concentrated on one cut [—Z2,Z2]- The transition 
from a two-cut to a one-cut regime is discussed in physical papers by Cicuta, 
Molinari, and Montaldi [CMM], Crnkovic and Moore [CM], Douglas, Seiberg, 
Shenker [DSS], Periwal and Shevitz [PeS], and others. 

Another proof of the limits (1.25) was recently given by Albeverio, Pastur, 
and Shcherbina (see [APS]) in a completely different approach based on some 
estimates of the Stieltjes transform of the spectral measure. 

A general ansatz on the structure of the semiclassical asymptotics of the 
functions ij^n{z) for a "generic" polynomial V{z) was proposed in the work [BZ] 
of Brezin and Zee. They considered n close to N , n = N + 0{1) , and suggested 
that for these ra's, 

(1.27) Vn(^) = cos(ArC(z) - {N - nMz) + x{z)) , 

for some functions f{z), C{z), ip{z), and x{^)- This fits in well the asymptotics 
(1.15), except for the factor (—1)" at x in (1-15), which is related to the two-cut 
structure of ipn{z)- 

Equation (1.7) also appears in the planar Feynman diagram expansions of 
Hermitian matrix models, which were introduced and studied in the classical 
papers [BIPZ], [BIZ], [IZ] by Brezin, Bessis, Itzykson, Parisi, and Zuber and 
in the well-known recent works by Brezin, Kazakov [BK], Douglas, Shenker 
[DS], and Gross, Migdal [GM] devoted to the matrix models for 2D quantum 
gravity (see also [Dem] and [Wit]). In fact, it is the latter context that broad- 
ened the interest to the Freud equation (1.7) and brought to the area new 
powerful analytic methods from the theory of integrable systems. It turns out 
[FIK1,2] that equation (1.7) admits a 2 x 2 matrix Lax pair representation 
(see equation (3.15) below), which allows one to identify the Freud equation 
(1.7) as a discrete Painleve I equation and imbeds it in the framework of the 
isomonodromy deformation method suggested in 1980 by Flaschka and Newell 
[FN] and by Jimbo, Miwa, and Ueno [JMU] (about analytical aspects of the 
method see, e.g., [IN] and [FI]). The relevant Riemann-Hilbert formalism for 
(1.7) was developed in [FIK1,2] as well. It was used in [FIKl-3] together with 
the isomonodromy method for the asymptotic analysis of the solution of (1.7), 
which is related to the double-scaling limit in the 2D quantum gravity studied 
in [BK], [DS], [GM]. 



194 



PAVEL BLEHER AND ALEXANDER ITS 



As a matter of fact, the solution of (1.7) which is analysed in [FIK] corre- 
sponds to the system of orthogonal polynomials on certain rays in the complex 
plane. Nevertheless, the basic elements of the Ricmann-Hilbert isomonodromy 
scheme suggested in [FIK] can be easily extended to an arbitrary system of 
orthogonal polynomials corresponding to a rational potential V{z) (cf. [FIK4]; 
see also Appendix E below). The proof of Theorem 1.1 is based on the ap- 
proach of [FIK] combined with the nonlinear steepest descent method proposed 
recently by Deift and Zhou [DZl] for analyzing the asymptotics of oscillatory 
matrix Riemann-Hilbert problems. We appeal to the Deift-Zhou method in 
Section 7 where we construct explicitly and then justify rigorously the as- 
ymptotic solution of the master Riemann-Hilbert problem associated to the 
orthogonal polynomials (see conditions (i)-(iii) and (5.16)-(5.18) below). 

The paper is organized as follows: 

In the next section we use the results of Theorem 1.1 for proving the 
universality of the local distribution of eigenvalues in the matrix model with 
quartic potential. We prove both the sinc-kernel universality at internal points 
and the Airy-kernel universality at the end-points of the support of the limiting 
spectral measure. It is important to note that for internal points the univer- 
sality was recently proved by a different technique in the paper by Pastur and 
Shcherbina [PS], for a general class of matrix models. 

In Sections 3-5 we reproduce in a slightly different way the results of 
[FIK] concerning the Lax pair representation of equation (1.7) and the matrix 
Riemann-Hilbert reformulation of the orthogonal polynomials. In particular, 
we show that there is an exact and simple relation between orthogonal polyno- 
mials Pn{z) and the 2x2 matrix- valued function ^n{z) which solves the fol- 
lowing matrix Riemann-Hilbert problem on a line (the problem (5.16)-(5.18) 
below) : 

(i) ^n{z) is analytic in C \ M, and it has a jump at the real line. 

(ii) *„(z) ~ (^^4) e-(^-"'""+^")'^^ z^oo, where 

An = ^ In /in, 0-3 = 

(iii) 4'„+(2) = *„_(2)S, Im3 = 0, S = 

As explained at the end of Section 5, in the setting of the Riemann-Hilbert 
problem (i)-(iii), the real quantities Rn and An are not the given data. They 
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are evaluated via the solution "^niz), which is determined by conditions (i)- 
(iii) uniquely without any prior specification of Rn and A„. Simultaneously, 
the function "ifniz) satisfies the Lax pair (see equation (3.15) below) whose 
second equation is the linear differential equation: 

(1.28) ^M^ = NAr,{z)^r,{z), 

where 

"^^^ V -I^\t + gz-" + g{Rn-i + Rn)] f + ^ + gzRn 

The jump matrix S in (iii) constitutes the only nontrivial Stokes matrix (for 
more details sec Sections 4, 5) corresponding to the system (1.28) with Rn 
generated by (1.4) and (1.5). This reduces the problem of the asymptotic 
analysis of the quantities Pn{z) and i?„ to the asymptotic solution of the 
matrix Riemann-Hilbert problem (i)-(iii), i.e. to the asymptotic solution of 
the corresponding inverse monodromy problem for differential equation (1.28). 

The short but important Section 6 provides a formal asymptotic ansatz 
indicated in (1.14) for the recurrence coefficients Rn. 

The proof of Theorem 1.1 is given in Section 7. The central point of 
the proof is a construction of an approximate solution ^^{z) to the Riemann- 
Hilbert problem (i)-(iii). To that end we consider an approximate differential 
equation (1.28) in which the coefficients i?„ are replaced by the numbers 

pO _ -t-{-ir^t^-A{n/N)g 
2g 

that are taken from our formal analysis in Section 6. The function *°(z) is 
constructed then as a semiclassical solution to this approximate differential 
equation with the large parameter A'^. Our semiclassical construction is based 
on the version of the complex WKB method which was recently suggested 
in [Kap] for asymptotic solution of the direct monodromy problems for the 
2x2 systems with rational coefficients. We partition the whole complex plane 
into several regions which separate different turning points of the approximate 
differential equation, and we construct WKB- and turning point semiclassical 
solutions in each of these regions (see details in Section 7). This provides us 
with an explicit matrix- valued function "^^{z) which solves asymptotically, as 
N — > oo, the basic Riemann-Hilbert problem (i)-(iii). Then we prove that the 
quotient '^niz) \^^{z)\ is equal to / + 0(A^^-'^(1 -|- Izl)"-*^) and this completes 
the proof of Theorem 1.1. We emphasize that we only use equation (1.28) 
to motivate our choice of the function ^'^(z). The uniform estimate for the 
difference, '^n{z) \^'^{z)\ ^ — 7, is proved by means independent of the WKB 
theory of differential equations. The main ideas and technique used in Section 
7 are based on the Deift-Zhou nonlinear steepest descent method [DZ]. 
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The Riemann-Hilbert reformulation of the orthogonal polynomials (1.1), 
(1.3) suggested in [FIK] plays the central role in our approach and in its ex- 
tention to the general rational potentials V{z). For that reason we decided to 
present with more detail the scheme of [FIK] in Appendix E. 

As mentioned above, the Freud equation (1.7) has a meaning similar to the 
discrete Painleve I equation. We refer the reader to the papers [FIZ], [NPCQ], 
[GRP], [Mag3,4], [Meh2] for more on the subject. As first noticed by Kitaev, 
equation (1.7) can also be interpreted as the Backlund-Schlezinger transform 
of the classical Painleve IV equation so that the coefficients i?„ coincide, in 
fact, with the special PIV function (sec [FIK1,3] for more details). This PIV 
function, in turn, can be expressed in terms of certain n x n determinants 
involving the parabolic cylinder functions (see [Mag4]). In this work however 
we do not use these algebraic connections to the modern Painleve theory. 
Instead, we use its analytical methods. 

The present paper is a revised and shortened version of our earlier 
preprint [BI]. 

2. Universality of the local distribution of eigenvalues 
in the matrix model 

Theorem 1.1 can be applied to proving the universality of the local distri- 
bution of eigenvalues in the matrix model with quartic potential. The matrix 
model is defined as follows. Let M = {Mjf;)j,k=i,...N be a Hermitian random 
matrix, with the probability distribution 

(2.1) i^NidM) = Z^^e-^T'-^(^)dM, 

where 

V{M) = ao + aiM + ■■■ + a2pM^P, asp > 0, 
is a polynomial, 

dM = Yl{d Re Mjk d Im M^-fc) JJ dMjj , 

j<k 3 

is the Lebesgue measure on the space of Hermitian matrices, and 

is the grand partition function. Let Ai < • • • < Ajv be eigenvalues of M. 
Consider the distribution function of the eigenvalues, 

FN{z) = N-^^#{j:\j<z}. 

and the density function 

Pn{z) = F'^{z). 
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In the matrix model we are interested in the following problems: 

(1) To calculate the limit density p{z) = liniN^ooPNiz)- 

(2) To calculate the limit local distribution (scaling limit) of eigenvalues at 
regular points, where p{z) is positive, and at end-points, where p{z) 
vanishes. 

(3) To calculate the free energy 

X logZN{ao,...,a2p) 
f{ao,...,a2p) = -hm^ 

and to find the points of nonanalyticity of / (critical points) in the space 
of the parameters oq, . . . , a2p- A further problem is to calculate the critical 
asymptotics of the recursive coefficients i?„ and of the local distribution 
of eigenvalues (double scaling limit) . 

Dyson [Dys] (see also [Mehl] and [TWl]) proved a formula which ex- 
presses the correlations between the eigenvalues of M in terms of orthogonal 
polynomials. Namely, the m-point correlation function is written as 

(2.2) KNm{zi, ...,Zm)= det {Qn^Zj, Zk))jf.^i^ „^ 

where 

N 

(2.3) Qn{z,w) = ^Vj(^)V'jM> 

i=i 

and ipj{z) is as defined in (1.9). When m = 1 the correlation function reduces 
to the function NpN^z); hence 

N 

p^{z) = N-'Y.^]{z). 

By the Christoffel-Darboux formula (see, e.g., [Sze]), the kernel Qn{z,w) can 
be written as 

^ / X V-^AT+l [tpN+l{z)ll^N{w) - lljNiz)lljN+l{w)] 

(2.4) Qn{z,W) = , 

z — w 

and 

.X ^/RN^[tp'N+li^)^Niz)-ll}'Niz)'tpN+l{z)] 

(2.5) pn{z) = — . 

The formula (1.24) is valid in a complex neighborhood of the interval 
[zi + 6,Z2 — 6] and this allows us to differentiate it. We will assume that 

t <tcT = -2-^5 
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(two-cut case); hence we can use n = N in the asymptotic formulae (1.15)- 
(1.18). For the sake of brevity we rewrite (1.24), (1.15) as 

Cz 



(2.6) 

where 
(2.7) 



i/jn = cos{N( + ri), 



C = C = C(^; A') = £ \Uoiv; \')\'/'dv + ^ ; 



v = -- 



dz 

(-ir 



(gz^ + tf 



x{z;X') = - 



arccos r, r = 



gz'^ + t 

2v^ ' 



and we drop terms of the order of N ^. In addition, (1.24) gives that modulo 
terms of the order of N~^, 



(2.8) 
where 
(2.9) 



daz;X') 



Cz 



cos{NC±^ -rj), 



1 



arccos q, 



gz'^ + t 



dX' 2 
The functions tpn satisfy the recursive equation 

Zl/Jn = \/Rn+l1pn+l + s/R^lpn-l 

(see (1.4)), hence from (2.6) and (2.8) wc obtain that 
zcos(A''^ — rj) cos(2?7) — zsin(A^C — i]) sin(2?7) 



2v^ 



= VRn+i cos{NC - rj) cos ^ - a/Rh+i sin(A^C - v) sin C 
+ ^/R^ cos{N( — rj) cos ^ + \/R^ sm{NC, — rj) sin ^. 

Equating the coefficients at cos(A/"C — rj) and sin(iV^ — ii)^ we obtain that 

(2.10) z cos 2r/ = Rn+i + /R^) cos 

zsin2r/ = {^J Rn+i - ^/Ru) svci^. 

These formulae can be checked directly from (1.14), (2.7) and (2.9). Differen- 
tiating (2.6) and (2.8) in z, we get that 



= -Czsin(iVC + 77) iV VCl + 0(1), 
<+i = -Czsin(iVC + e - »?) iV^G + 0(1); 
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hence by (2.5), modulo terms of the order of A'^ ^, 
(2.11) 



PN = ^fR^i Cz" [- sm(iVC + C - r/) cos(iVC + r?) 

+ cos(ArC + ^ - ry) sin(A^C + 
= \/ Rn+1 C^z^ sin(277 — ^) = ^/Rn^ C^z^ {sin 2r] cos ^ — cos 27] sin 

and by (2.10), 

PN = \/ Rn+1 C^z ( Rn+1 - \/Rn) sin ^ cos ^ 
-(-^i^AT+i + ^/Rn) sin^cos^ 
= -^/Rn+iRn C^z sin2^. 
Since modulo terms of the order of A^~^, 



Rn+iRn = - ; 

9 

we obtain that 



IT 



sin 2^ = sin(— arccosg) = — -\/l — q^, 



V9 



PN = ^zVi~q'^ + o{N-^] 

Substitution of the value of q gives that 

PN{z)=piz) + OiN-^), 



where 
(2.12) 

p(z) = 3l?7o(z;l)|V2 

and 
(2.13) 



TT 



\Z\ 
TT 



1/2 



9\z\ 
27r 



^1,2 



9 



This gives an explicit formula for the limiting density p = p{z) of eigenvalues 
(integrated density of states). In a completely different approach, based on the 
Coulomb gas representation of the matrix model, this formula is derived in the 
work [BPS] of Boutet de Monvel, Pastur, and Shcherbina, as an application of 
the proven (in [BPS]) variational principle for the integrated density of states. 

The scaling limit of the correlation function KNmizi, ■ ■ ■ ,Zm) Sit sl regular 
point z, where p{z) > 0, is defined as 



Kmiui,..., Um) = lim [Np{z)] '^Kntu ( z + 
N—>-oo \ 



Ul 

Np{z) 



,...,z + 



Np{z) 
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Observe that Km{ui, . . . ,Um) is the hmiting m-point correlation function of 
the rescaled eigenvalues 

fij = Np{z){Xj - z). 

The rescaling reduces the mean value of the spacing /Xj+i — nj to 1. Prom 
Dyson's formula (2.2), 

(2.14) Km{ui, . . .,Um) = det{Q{uj,Uk)).j^^-^^ 

where 



(2.15) Qiu,v)= lim [Np{z)] ^Qn(z + 

AT— »oo y 

By (2.4), 

(2.16) [jVpW]-Q„(^. + + _ 



U V 



Np{z)' Np{z] 



■Tn \ z+ ———-,z + 



u — V \ Np{z) ' Np{z) J ' 
where 

Tn{z,w) = iPn+i{z)'>Pm{w) - 1pN{z)'lpN+l{w)- 

By (2.6) and (2.8), modulo terms of the order of N~^, 
(2.17) 

, / u \ Cz . CzU 

tpN \ z+ = —7= cos[NQ + a + r]), a = 



hence 
(2.18) 



Tn[z+—^,z + 



Np{z) ' Np{z) 



[cos(ArC + a + C - r?) cos(ArC + /3 + ??) 

Cz 

- cos(A^C + a + ?y) cos(A^C + /9 + C - ??)] 
C'^z^ 

[cos(a + ^ - /3 - 277) - cos(a - ^ - /3 + 2r?)] 



sin(2?7 — ^) sin(a — /3), 



where 



r9 1Q^ l^0(^)l^/'^ n 
(2.19) «=—— = - , , ,^ = TTU, p = TTV . 
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By (2.11) and (2.12), 

V^^CV sin(2r? - = p{z) = -y/\Uoiz-l)\ = U^; 1); 

TT 

hence (2.18) implies that 

V^n+iTn + J^p(^'^+ J^p(^] = V^N+i sin(2r? - ,e) sin(a - /?) 

sin(Q; — /3) sin 7r(n — v) 



TT TT 



and, by (2.15), (2.16), 

Q{u,v) 



sin7r(M 



^{u — v) 



This proves the Dyson sine-kernel for the local distribution of eigenvalues at a 
regular point z. In a completely different approach, the sine-kernel at regular 
points is proved in [PS]. 

Remark. It follows from the Dyson sine-kernel, due to the Gaudin formula 
(see, e.g., [Mehl]), that the spacing distribution of eigenvalues is determined 
by the Fredholm determinant det(l — Q{x,y))x,yeJ- The asymptotics of this 
determinant as | J| — oo has been studied intensively since the classical works 
by des Cloizeaux, Dyson, Gaudin, Mehta, and Widom (see [Mehl] for the 
history of the subject). The Riemann-Hilbert approach to this asymptotics 
has been developed in the paper [DIZ]. 

At the endpoints of the spectrum we use the semiclassical asymptotics 
(1.18), and it leads to the Airy kernel (cf. the papers of Bowick and Brezin 
[BB], Forrester [For], Moore [Mo], and Tracy and Widom [TW2], where the 
Airy kernel is discussed for the Gaussian matrix model and some other related 
models, and, in addition, some nonrigorous arguments are given for general 
matrix models). Consider for the sake of definiteness z = Z2- 

By (1.18), 

(2.20) Vn = ^^ M{N^'^w)+0{N-^) , D = ^, 



where w is defined as in (1.19). From (1.19), 
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This allows us to derive from (2.20) that 
(2.21) 



1pn±l 



Ai (7V2/3^Q + N-^l^u) + 0{N-^) 



Ai (A^2/3^^ ^ iV-V3p _ iv-V3^^ + 0{N-^] 



where 



(2.22) '^0 = ^Pq{z;>!) 

p = p{z; A') = 



3 

e(z;A' 



and 
(2.23) 



e(^;A') = - 
^(^;A') = - 



cosh^^ g 



, \ 2/3 

y'Uo{v;\')dv\ , 
u = u{z; A') 



2 

(-ir 



9 



+ t 



cosh ^ r, r 



2y/Yg-tq 



4 2VA^g - 1 

The formulae (2.22), (2.23) define the functions Lpo{z;\'), p{z;\') and uj{z; 
for z > Z2- It is easy to check that these functions are analytic in ^; at ^; = 
and they can be continued analytically to the interval z > zi. In addition, 



(2.24) Uo{z2;\') = 0, 



dUo 
dz 

dz 



(.2;A') = x = 2(A')^/V/%^ 
(.2;A') = xV3 = 2V3(A')VV/^z2; 



p{z2) = -2-V3(A')-V3. 

a;(.2) = -^2V3(A')-V3,,,-i. 

We will consider 

(2.26) z = Z2 + N-^/^a, w = Z2 + N'^/^P, 

where a and (3 are fixed. 

Substitution of (2.21) into the recursive equation 

Zipn = V -Rn+lV'n+1 + \/^'0n-l 

gives the equations 

(2.27) Z2 = ^/R^l + ^/R^, 



Z2UJ = ^/Rn+l{p -U}) + V^(-P - 1^), 
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from whence 

(2.28) {y/Rn+i + VRu) 2w = ( V^n+i - VRu) P- 

Similarly, 



hence 

(2.29) 2u = tllllP , 

which agrees with (2.24). 

Substituting the formulae (2.21) into (2.4) and throwing away terms of 
the lower order, we obtain that 

(2.30) 



Qn{z,w) 



{z - w)(p'o 



1/3 



Ai {N^/^^oiz) + N-^l^p - N-^/^uj) Ai (iV2/Vo(«^) + iV" 
- Ai (Ar2/3^o(-2) + N-^/^u) Ai (Ar2/3<^o(w^) + N-'^'^p - N'^/^u) 

where <^q, p and u are taken at Z2- Taking the linear part of Ai we obtain 
that 

QNiz,w) = [Ai(«)Ai'(^) - Ai'(u)Ai(^)] (2a;-p)iV-V3, 

{z - w)(p'o 

where 

u = ifQa, V = ipQ 13. 
By (2.28) and (2.24), modulo terms of the order of N~'^/^, 

hence 

Q^(., ^) = ArV32V3,i/2,^ Ai(.)Ai»-Ai-(.)Ai(.) ^ ^ /3 

u — 

Thus, 



1 / u \ _ Ai(u) Ai'(i;) - Ai'(?x) Ai(t;) 



lim 
where 



c = ^(,(z2;l) = 21/351/2^2. 



This proves the Airy kernel at the endpoint Z2- The endpoint zi is treated 
similarly. 
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3. The Lax pair for the Freud equation 



Let 



(3.1) 

Then 

(3.2) 



■■Pn{z)e 



-NV{z)/2 



/oo 
-oo 



A recursive equation for tpniz) follows from (1.4): 



(3.3) 

In addition, 
(3.4) 

V'U^) 



+ 
+ 



Zl/Jniz) = Rl/^itpn+l{z) + Rl/'^tpn-l{z). 

N^Rl/^,Rl%Rl%)^n+3{z) 

N I + iV I Rl!+l{Rn + Rn+1 + Rn+2) 

N^-Rl/^ + N^ Rl!^{Rn-l + Rn + Rn+l) 



i^n+l{z) 
i^n-liz) 



[NlRT-.RT-.Rl!')i^n-.{z). 



Let 
(3.5) 



^^niz) 



i^n{z) 
V'n-l 



Then combining (3.3) with (3.4), one can obtain (cf (3.1-7) in [FIK2]) that 

J ^n+l{z) = Un{z)^n{z), 

\ ^'^{z)=NAn{z)^^{z), 



(3.6) 

where 
(3.7) 

and 
(3.8) 

An{z) = 

Observe that 



Un{z) 



-"■n+1 Z 



R!rl'^[t + gz'^ + g{Rn + Rn+l)]\ 
-R^'' [t + + g(^R^_^ + R^y^ f + ¥ + 9zRn ) 



[i^ + a^+gzRn, 



tr^„(2) =0 
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and 
(3.9) 

deiAn{z) = - I TT + -?r +9Rn (t + gRn-l + gRn + gRn+l) + RnOn-lOn-, 



where 

(3.10) en = t + gRn+gRn+i. 
Due to (1.7), we can rewrite detAn{z) as 

(tz z'^\'^ nz'^ 
y + ^j +^+RnOn-lOn- 

The compatibiUty condition of equations (3.6) is 

(3.11) Kiz) = NAn+l{z)Un{z) - NUn{z)An{z). 

Restricting this equation to the matrix element U!^ u{z) we obtain that 

(3.12) Jn+i — J„ = 1, 
where 

Jn = NRn[t + g{Rn-l + Rn + Rn+l)]- 

Hence Jn = n + const. Since Jq = 0, in fact, J„ = n. This means, that the 
compatibihty condition (3.10), together with the initial value Jq = imply 
(1.7), and thus equations (3.6) give the Lax pair for the nonlinear difference 
Preud equation (1.7). In addition, equation (3.12) gives the recursive equation 

(3.13) Rn+l9ndn+l = R-ndn-ldn + • 

Equations (3.6) have two linear independent solutions and '^n{z) is one of 
them. We will consider another solution, 

^n{z) 
^n-l{z) 



and the 2x2 matrix 



(3.14) ^.(.) = f/"(;\ 

which satisfies the same equations, 

f "^n+liz) = Un{z)'i>n{z), 



(3.15) 



<(z) = NA^{z)^^{z). 



To define $n(z) we consider an arbitrary, linearly independent with ^i(z), 
solution of the differential equation ^[{z) = NAi{z)^i{z), and then define 
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^^(z), n > 2, with the help of the recursive equation ^n+iiz) = Un{z)^niz)- 
Equation (3.11) then leads to the differential equation ^n{z) = NAn^n{z) 
for n > 2. Note that the equation ^n+i{z) = Un{z)^n{z) means that ipn{z) 
satisfies the recursive equation (3.3); i.e., 

(3.16) Zipniz) = R%^n+l{z) + R]i\n-l{z). 

Since tr ^^(z) = 0, the second equation in (3.15) implies that 

(3.17) det^„(z) = C/0 
is independent of z; i.e., 

(3.18) ^n{z)^n-l{z) - ^n-l{z)^n{z) = C, C = C{n). 

This enables us to derive a first order differential equation for ipn{z). Namely, 
from the first row in '^n{z) = NAn'^n{z) we have that 

(3.19) ip'^{z) = Nan{z)ipniz) + Nai2iz)ipn-iiz), 
<(z) = Nau{z)^n{z) + Nai2{z)ilJn-i{^), 

where An{z) = {aij{z)) ■ hence 

'^n{z)(p'n{z) - (pn{z)^Pn{z) = Nai2{z)[ll)n{z)ipn-l{z) - V'n-1 (^)'/'n (2)] 

= CNai2{z), 

and 

^3 20) ,^ipniz)\' _CNai2iz) 



i^n{z)J i^liz) 

In a similar way from the second row in (3.15), ^'n{z) = N An{z)'^ n{z) , we 
obtain that 

,3.21) (^y._^^. „>i. 

It is useful to note that (3.18) allows us to express ipn-i{z) in terms of ^Pn{z): 
(3.22) <^^_^(^) = ^!lziM^4^)+ ^ 



•^^niz) 1pn{z) 

Remark. The system of two differential equations of the first order, 
= NAn'^n, can be reduced to one equation of the second order (cf. [Sho]). 
Namely, from the first equation of the system we can express V'n-i in terms 
of ipn, 

(3.23) = iv-i— <-^Vn, 

ai2 ai2 
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and then we can substitute this expression into the second equation of the 
system, which gives 



(3.24) - ^ + N\aua22 - 012021) Vn - iVai2 

ai2 

With the help of the substitution 

(3.25) Vn = al^^Cn 
we reduce (3.24) to the Schrodinger equation 

(3.26) -C + N'^UCn = 0, 

where 
(3.27) 

U = -((111022 - 012021) + I a[i - an — 



Oil 

012 



012 



N- 



ii2 3(a'^2)^ 



2ai2 4af2 



By (3.8), 
(3.28) 



tz gz 

-On = 022 = Y + — + gzRn, 



au = Rl/H0n + gz"), 021 = -Rl/\en-i + gz 

which gives 
(3.29) U{z) 



g^z^ tgz'^ 



f n 



■N- 



't ^ 3gz^ gz^{t + gz^ + 2gR,,) 

^ + ^+gRn ^^^^ 



+ N- 



g(2gZ^ - 0n) 

. igz^ + 9nr 

It is convenient to write U {z) as 



(3.30) 

where 
(3.31) 

Uo{z) = z^ 



Uiz) = Uoiz) + Ui{z) + U2iz), 



TV ' 



Uiiz) = N-' (^l + gRn^ , 

U2iz) = -Rnen-ien-N-^ 



0n{t + gz'^ + 2gR, 
gz^ + On 



+ N- 



' g{2gz^ - 9, 
[{gz^ + Ony 
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4. The Stokes phenomenon 

Consider the sectors 

(4.1) 

nj = lz:-+ ^ > +e<Bxgz<—+ ^ ^ -gj, J = 1,2,3,4, 
g: > 0, on a complex plane where the function 

(4.2) ^„(z) = ^P„(z)e-^(l^'+^') 
goes to infinity as z — oo. Let hj be the bisector of Oj, 

(4.3) bj = {z: z = rujj, r>{)] , = e'(^+^^^) , j = 1,2,3,4. 

For a fixed n > 1, let (7„ > be a constant which will be chosen later, and let 
us consider the following special solution to equation (3.20) with C = Cn- 

(4.4) 

^„,.(.) = CnNMz) r ^ dn = CnNhl/'Pniz) e-^(^'+i^') 



Then ipnj{z) — as z — oo in Evaluating the integral in (4.4) with the 
help of the Laplace method, we obtain the asymptotic expansion of <Pnj{z) in 



(4.5) 

/ oo \ 



z^ 

k=l 



Let us define the function (pn^ij{z) by the formula (3.22), 

(4-6) (^„_ij(z) = <yfni(^;) + 



Then the vector-function 

"'^ ^ V'^n-ij(2)y 

satisfies the differential equation 

^U^) = NA^{z)^nj{^). 

Hence the function ipn-ij{z) satisfies the first-order differential equation (3.21) 
with C = Cn- In addition, (4.6) gives that <^„_ij(z) — > as 2; — > oo in ilj. 
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This implies that 

a2i{u) 



(4.7) (Pn-l,jiz)= - CnNlPri-liz) / 

J to 



.,,2 

WoOO 



du 



I 



UJjOO 

which for n > 2 can be rewritten as 

,1/2 



r Rn-At + 9U' + g{Rn-l + Rn)] ^NiW+\u^)^^ 



The last expression is similar to (4.4) if we take 

Cn-l = CnRn'^Rn-l : n > 2, 

or 

Cu-iRH-i = CnRn"^, n> 2. 

We solve the last equation recursively. The initial constant Ci is a free param- 
eter. We put Ci = R^^^"^. This gives Cn = Rn^'"^ , so that from (4.4) and (4.7) 
we derive 

(4.8) ^„,(.) = iV/.y2p„(,)e-^(K+f^^) 

^ r t + gu^+g(i^ +i?„+,) ^^(|„^+K)^^^ 

and 

(4.9) ^„.(,,.ftyV»-.e^('i + f;2|V 
Observe that 5nj(-2) satisfies the recursive equation 

^n+lA^) = Un{z)^nj{z). 

Indeed, consider the vector-function 



Then, since ^>„j(z) satisfies the differential equation 



fnj{z) 
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the equation (3.11) imphes that ^^_^_ij{z) satisfies the differential equation 

the same as for ^n^ij{z). In addition, ^(z) ^ as 2; — oo in J7 j . Hence 
^n+i,j(^) coincides with $^+1^(2;) up to a constant factor, but their second 
components are the same; hence ^(z) and ^n+i,j{z) coincide, which was 
stated. 

It is also worth noticing that ^nji^) as a solution of linear ODE with 
polynomial coefficients is an entire function of z. 

As a matter of fact, the Laplace method combined with a deformation 
of the contour of integration gives the asymptotic expansion (4.9) in a bigger 
domain, 
(4.10) 

„ (it ttH — 1) Svr nii — 1) 1 

Si = |-g+ '\ ^ +£<argz<— + ^ -ej, e > 0. 

Thus, the matrix- valued function 

(4.11) *.«=(,!:';;!, ^}^;) ^ 

is an entire function of z, and according to (4.2) and (4.9) it has the asymptotic 
expansion 

(4.12) ^^j{z) ~ ^) ^-(-^-nln.+A„).3^ ^ ^ ^ ^ 5..^ 



1 
-1 



where 

(4.13) (73 = 

is the Pauli matrix, 

(4.14) A„ = ^ln/i„, 

and Ffe are some 2x2 matrices which only depend on n. Prom (4.2) and (4.9) 
we obtain that 

The independence on j of the coefficients in the asymptotic series 
(4.12) constitutes the so-called Stokes phenomenon, and the functions ^nj{z) 
form the set of canonical solutions (see, e.g., [Sib]) for the equation '^'{z) = 
NAn{z)^{z). The corresponding Stokes matrices (see equation (5.3) below) 
are evaluated in the next section. 
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5. The Riemann-Hilbert problem 

Since both ^nj{z) for different j and ^n(^) satisfy the same differential 
equation (3.6), they are linearly dependent, 

where j is defined mod 4. The domains Hj and Sj+i intersect and in the 

intersection the functions ^nj{z) and ^nj^i(z) grow to infinity and have the 
same asymptotic expansion (see (4.12)). On the other hand ^^(-z) goes to zero 
in this intersection, hence q = 1, so that 

(5.1) ^n,j+liz) = ^nj{z) + S*n(z). 

Since both ^njiz) and 'i'niz) satisfy the same recursive equation (3.6), the 
coefficient s does not depend on n, but in general it depends on j, s = sj. The 
equation (5.1) implies that 

(5.2) ipn,j+l{z) = ^Pnj{z) + Sjll^niz). 

We can rewrite (5.1) in matrix form as 

(5.3) ^n,j+l{z)=^nj{z)Sj, 

where 

(5.4) S,= (J ^1). 

To determine Sj consider (5.2) at n = 0. Prom (4.8) it follows that 

(5.5) ^ojiz) = iV/iJ/'e-^d^'+i^') {t + gu" + 5i?i)e^(^'+f"')dn. 
Putting z = in (5.2) we get 

_ yo,j+i(0) - </^0i(0) 
^o(O) 

Since 

V'o(O) = h-^'^ 

(see (4.2)) and 

¥Poi(0) = Nhl'^ f {t + gu^ + 3i?i)e^(t«'+f«') du 

J UJjOO 

we obtain that 

(5.6) sj = Nho / {t + gu^ + gRi)e^\^'' +f " ) du, j = 1, 2, 3, 4. 
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The change of variable u ^ —u gives 

(5.7) S3 = -Sl, S4 = -S2- 

Another way to compute Sj is to use the Cauchy type integral. 
The function 

_ NV(z) 

(5.8) y,y(z)=e- 2 99„j(z), i = 1,2,3,4, 
is an entire function of z and in it has the asymptotics 

(oo 

In addition, 



l2k 
2k 



-NV(z) 

(5.10) yn,j+i{z) = ynj{z) + sje 2 ip^iz). 

Observe that the domain Hj contains the j*^ quadrant, 

A, = |^:^— <arg.<-j; 

hence (5.9) holds in Aj. This allows us to solve (5.10) with the help of the 
Cauchy type integral. Namely, 

(5.11) 

NV(u) ■ NV(u) 

roo g ^^„(«) ^ Sl e —Mu) 



(z) =— n ^ ' ^"^^-^ du+— r 



27ri ./-oo u — z 27ri ./_,-,^ « — 2; 



where yn(^) is a piecewise analytic function which coincides with ynj{z) in the 
quadrant Aj. Expanding 



It — Z 2 ^ — ' 

k=0 



we obtain from (5.11) the asymptotic expansion of yniz) as z ^ 00: 

(5.12) yn(^)~-2^ -k , 



z 

k=0 



where 



1 /"°° 

(5.13) ank = -TT^ / /i-^/'Pn(«)u'e-^^W du, 

bnk = / /i-V^Pn(u)u^e-^^W du. 

2m J-ioo 
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Since by (5.9), y„(z) = 0(z-"-i), 

SAttnk + Sibnk = 0, = 0, 1, . . . , n - 1, 

but 

O'nk — foi" k — 0, 1, . . . , n — 1 in virtue of the orthogonality property (1-3); 
hence sib^k = for these A;'s. Let us take n = 2 and A; = 0. In this case 

/oo 
-oo 

is obviously positive; hence 620 7^ 0. This implies 
(5.14) si = S3 = 0. 

By (5.2) this means that 

Let us find S4. From (5.9) we know that 

y^{z) = hy\z-' + ...)■ 

1/2 

hence by (5.12), ^4000 = • -'■^ addition, by (5.13), 

1 1.1/2 

aoo = -^/.o • 

This gives 

S4 = —2-Ki. 

It is interesting to notice that the values of si = S3 = and S2 = S4 = —2m 
can be obtained directly from formula (5.6) as well, although the calculations 
are more involved. Namely, the recurrent coefficient R\ is expressed in terms 
of the parabolic cylinder function: 

1/2 



z=t. / ^ 



Then, using two different integral representations of the function D_ii2{z) (see, 
e.g., [BE]), 



poo 



\^ Jo 
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one derives from (5.6) that si and S4 are expressed in terms of the Wronskians, 



Si = V^TT 

S4 = —V^iri 



0, 



-2m, 



as stated above. 

Now we can formulate the Riemann-Hilbert problem. Define 
(5.15) 

^'n-(^) = 



Let 



'^kn{z) <Pnl{z) 
'(pn-l{z) ipn-l,l{z) 



1pn-l{z) (Pn-l,3{z) 



^niz) 



^„+(z), if lmz>0, 
^n-{z), if lmz<0. 
Then '^n{z) has the asymptotic expansion 



(5.16) 



^k=0 



z 00, 



where A„ = ^ In /i„ and 



(5.17) To 

On the real line 

(5.18) 

where 

(5.19) 

In addition, 
(5.20) 



R 



-1/2 , 







^n+iz) = 'i>n-{z)S, lmz = 0, 



1 -2m 
1 



det^n{z) = R-^/^ ^0. 



It must be emphasized that in the setting of the Riemann-Hilbert problem 
(5.16)-(5.18) the real quantities i?„ and A„ are not the given data. They 
are evaluated via the solution "ifniz)- Indeed, suppose that "ifniz) is another 
function satisfying (5.16)-(5.18) with perhaps another Rn and another A„. 
Consider the matrix ratio: 

X*{z) = [e^-''^r^'^n{z)][e^-'^^r^'^n{z)]-\ 
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Since the jump matrix S for both ^'n(-z) and ^n(-?) is the same, function X*(z) 
has no jump across the real hne. Therefore, it is an entire function equal to / 
at 2; = oo. Hence 

X*{z)^I, 

or 

Substituting into this identity asymptotic expansion (5.16) and equating the 
terms of order we come up with the matrix equation, 

or 

/ e2^« \ _ f e^^" \ 
[Rne-^^" J ~ {Rne-^'^" J ' 

hence 

Aji = and Rn = Rn- 

The Riemann-Hilbert problem (5.16)-(5.18) for the orthogonal polyno- 
mials (1.1), (1.3) first appeared in [FIK2,3] as a particular case of the more 
general Riemann-Hilbert problem related to the system of orthogonal polyno- 
mials on the cross R|JiIR with the same weight function . Also, in [FIK2,3] 
an alternative approach to the derivation of the basic equations (5.15), which 
can be easily extended on an arbitrary system of orthogonal polynomials (cf. 
[FIK4]), was suggested. We present the method of [FIK] in Appendix E. 



6. Formal asymptotic expansion for Rn 

To guess a behavior of Rn let us evaluate Rn for small n. We have i?o = 

and 

ho /-e-^(^^-^^)^dz " 
Using the change of variable 



we obtain that 



jryjl + any/'e-^^du (2^)1/2 

ill = 5 7 = T', 
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Expanding (1 + au)^^^"^ into Taylor's series and evaluating the Gaussian inte- 
grals we derive that 

3a^\ \t\ 1 3g 



Ri = MUl 



a 

T 



+ 0(N-'') = ^ - 



g \t\N |t|37V2 



+ 0{N-''). 



Now from the Freud equation (1.7) we obtain subsequently that 



R2 



+ 



4<7 



+ 0(iV-3), 



t\N \t\^N^ 



R4 = 



\t\N 



+ 0(iV-2), 



\t\N 



+ 0{N- 



These formulae suggest that for a fixed n as iV — 00, 



n 



(\t[ 

9 m 



+ 0{N- 



if n = 2k + 1, 



n 
JtlN 



+ 0{N- 



if n = 2k. 



n 



More generally we expect that if — is bounded and N ^ 00 then 



(6.1) 



' (-^) + 0(Ar-2) n = 2k + l, 
L (^^'^ + 0{N-^) if n = 2k, 



where i?(A) and L{X) are smooth functions. Substituting this ansatz into the 
Freud equation (1.7) and neglecting terms of the order of 0{N~^), we obtain 
the equations 

(6.2) X = R[t + g{2L + R)], 

X = L[t + g{2R + L)]. 

Equating the expressions on the right we obtain that 

(R-L) [t + g{R + L)] = 0. 

We assume that R L; hence 

(6.3) t + g{R + L) = 0. 
Combining this with (6.2) we obtain 

(6.4) A = gRL, 

so that R, L are solutions of the quadratic equation 



(6.5) 



+ - u + - = 0, 



9 



9 
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which are 

(6.6) i?,L= . 

Actually, we can find a formal asymptotic expansion of Rn in powers of 
Put 

{ R{n/N) if n = 2fc + l, 
^ ~ j L{n/N) if n = 2k. 

Then (1.7) is equivalent to 

(6.7) X = R{t + gR + 2gL + N-^gAL), 
X = L{t + gL + 2gR + N-'^gAR), 

where 

.f... _ /(A-^)-2/(A) + /(A + ^) 

Let us substitute the expansions 

L(A) = Lo(A) + N-^Li{X) + ...; R{X) = i?o(A) + N'^RiiX) + ... 

into (6.7) and equate coefficients at powers of N'^. Then the equations on 
Lo(A), -Ro(A) coincide with (6.2), hence Lq{X), Ro{X) are given by (6.6). Equat- 
ing coefficients at we obtain the system of equations 

{(i?o + Lq)Ri + 2RqLi = —RqALq, 
2LoRi + {Ro + Lo)Li = -LqARq. 
Solving this system we find Ri,Li and so on. 

7. Proof of the main theorem: 
The asymptotic Riemann-Hilbert problem 

In this section we construct explicitly, using semiclassical formulae, a 
piecewise analytic function ^^(z) which is an approximate solution to the basic 
Riemann-Hilbert problem (5.16)-(5.18). Then we prove that "^^(z) is close to 
*n(^) (the Deift-Zhou method, see [DZ]). 

Define the numbers (cf. (6.6)) 

(7U po_ -t-(-l)"(t^-4Aff)V^ n 

(7.1) K- , A--, 

_t+{-ir{t'-4Xg)y' 
^"±1 " 2~g ' 

and consider an auxiliary matrix differential equation, 

(7.2) *'(z) = NA\z)'^{z), 



(6.8) 
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where the matrix 

is defined as 
(7.3) 



_ (-l)"(t2_4Aff)V2 



(tz gz^ o\ gz^ 

an = - ( y + -Y+3zRn] = anz- — , 

a?2 = {<f"9z\ 

The matrix ^^(2;) is obtained from the matrix An{z) by the substitution of i?° 
for Rj, j = n—1, n, n+l. The function ^^{z) is constructed as a semiclassical 
approximation to the differential equation (7.2). 
A direct calculation gives that 

- del A\z) = ^ [it + 9z^f-4Xg]. 

Define the //-function as 
(7.4) 

/x(z) = V-detAO(z) = ^it + gz^y^AXg = (z^ - zl){z^ - z^), 

where we take the branch of the square root which is positive on the positive 
half-axis. The zeros of the //-function are the turning points for the equation 
(7.2). There are five turning points: a double turning point at 2; = and four 
simple turning points at ±zi and ±2:2 where 



^1,2 

We assume that 



e < A < e 

4<7 



(cf. (1.11)), which ensures that both zi and Z2 are real and separated from 
zero. 

Let n be a neighborhood of the interval [21 , 2:2] . For the sake of definiteness 
we define O as an ellipse with foci at zi and 2:2, 

f (x — 20)^ 

\l= < z = X + iy: 7^ \- tk < I 

Zl+ Z2 2 / \2 , i2 

20 = — ^ — , a = [zq - zi) +b . 
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We will assume that 6 > is independent of N and it is sufficiently small so 
that Q, does not contain the origin. Let 

r^i = n {Rez < zo}, Q2 = ^^{Rez>zo}, 

n"" = nn{imz>o}, ff^ = nr\{imz<o}, 

and assume that the boundaries dQ,, dQ,i, and dQ^ are positively oriented (see 
Fig. 2). Denote 

Ci = dnn {Re z < zo}, C2 = dnn {Re z > zo}. 

Let vi^2 = zo±ib be two vertices of the ellipse Q and 

lo = [vi,V2], lo=[vi,Zo], ^0 = [2:0,^2]. 

Wc shall also take the following convention: given any set A C C we denote by 
A" and A'^ the closed upper and lower parts of set A with respect to the real 
axis, so that 

J7« 2 = ^1,2 n {Im z > 0}, nf 2 = 5^1,2 n {Im z < 0}, 
Ci,2 = Ci,2 n {Im z > 0}, Ct^2 = Ci,2 n {Im z < 0}, 

etc. (see Fig. 2). 

Our plan is the following: 

(1) Outside of the domain fiU (— where — f2 = {z : — z G fi}, we define 
^°(z) by a WKB formula. 

Im z 




Rez 



Figure 2. The set £2 
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(2) In the domains O^'g we define ^^(z) by turning point formulae, with 
the help of the Airy function. In we define ^"(z) = (-l)"(73*°(-^;)(73. 

(3) We check that the multiplicative jump matrix of '^^{z) between 
and coincides with the matrix S in (5.19), and the multiplicative jump 
matrix on di^ is close to the unit matrix, as well as the multiplicative jump 
matrix between and ^72- 

This will give us an approximate explicit solution ^'^(z) to the Riemann- 
Hilbert problem (5.16)-(5.18), and then we will establish the closeness of ^^{z) 
to *„(z). 

As the first step of our plan, we define ^'°(z) outside of ^7 by a WKB 
formula. To get this formula we diagonalize iteratively the equation (7.2). 
The diagonalization is described in detail by Kapaev [Kap]. We will use the 
WKB formula with the error term 0{N''^). In this approximation we get that 
(7.5) ^^{z) = *wkb(^) = Cor(z)e-^^W'^3-r(z)+Cia3^ 

where 



(7.6) T{z) = "-^"^ , C{z) = £ t^iu) du, 

t{z)= ! di&.gT-^{z)T'{z)dz, 

J zo 



and 



(7.7) M^) = V-det^O(z)= {a\^{z)y ^a\^{z)a\^{z) 



1/2 



In (7.6) diagM stands for a diagonal matrix which coincides with M on the 
diagonal. The formula (7.5) gives an approximate solution of the differential 
equation (7.2), up to the terms of the order of iV~^, and this is a motivation 
for us to define '^^{z) outside of J7 by (7.5). However, it is important to 
underscore here that our ultimate goal is a construction of an approximate 
solution to the Riemann-Hilbert problem (5.16)~(5.18) and we are not going 
to prove any results concerning the differential equation (7.2). It serves only 
as an auxiliary tool in our construction of '^^{z). 
Observe that 

IJ,{z) = z^/XgsJx^ - 1, X = ^J" £^ , 

2y/\g 

and the values z = zi and z = Z2 correspond to x = — 1 and x = 1, respectively. 
We fix the branches of the multivalued functions in (7.6) by the inequalities 



— — < argz < — , \z\ > Z2, — TT < arg va; — 1 < tt. 



TT < arg(a; + \/ x'^ — 1) < tt. 
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This defines n{z) as an analytic function on C, witli two cuts [—Z2, —zi] and 
[zi,Z2]- The function ijl{z), in addition, satifies the symmetry equations, 



fi{-z) = -fi{z), n{z)=fi{z). 

The formula for fj,{z) leads to the following expression for the function ^{z): 
(7.70 

A r n; 7 , . r-^ 7.1 _ * + 9^"^ _ 



^(z) = — x\/x^—l — ln(x + \/ x'^ — 1 



The function ^{z) is analytic and it has one cut (— oo, Z2]- 
The matrix T{z) diagonalizes the matrix A^{z), 

T-\z)A'{z)T{z) = -fi{z)as, 

so that TiJ'iz) are eigenvalues of the matrix A'^{z) and T{z) is the matrix 
of corresponding eigenvectors. The matrix diagT^^(z)T'(z) turns out to be 
proportional to the unit matrix and a straightforward calculation gives that 
up to a nonessential constant factor. 



2/i y ' \o 1 

(see Appendix A below). Thus (7.5) reduces to 

(7.9) *°(z) = ^'wkb(^) = Co^o(z)e-^«(^)-3+^l-^ 

where Co 7^ 0, Ci are constants and 



(7.10) To{z) = 



u „0 X 1/2 / 1 



,0 
-'12 



2/x 

We shall fix the branch of the square root by the condition, 

> 0, Z > Z2. 



\ 1/2 



2fi 

It is easy to check that Tq{z) has no singularities except for those at 2; = ±zi 2, 
and that 

detro(2) = 1. 
The last equation in turn implies that 

det*wKB(^) =C^ 7^0. 

Consider the region 

c* = c\{nu{-n)). 
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oo. 



Lemma 7.1. The function ^'wkb(-2) defined in (7.9), is a single-valued 
analytic function on C* , which, in addition, satisfies the symmetry equation 

(7.11) *WKb(^) = (-l)V3*WKB(-^)tT3. 

Proof of Lemma 7.1 is given in Appendix B below. As z — > oo, 

(7.12) T„(.)=(j ;)+^-'((^';v. ^'^"^od^i-). 

In addition, from (7.6), (7.7'), 

(7.13) 5(.) = IM-!l]^+^„ + 0(N-), , 
where 

t^ X g X 

(7.14) = 
Thus, 

(7.15) *0(z) ~ (f^ ^] Coe-(^-'^i°^+^^"-^0-3, 



oo, 



^k=0 



with 

(7.16) «o-(^o ij, Oi- 1^^^0)1/2 

Turning Point Functions. To construct ^'^{z) in the sets Oi_2 we use 
turning point functions which are expressed in terms of the Airy function. Let 
us consider first the set ^2- On this set we are looking for a solution to (7.2) 
in the form 

(7.17) "^{z) = W{z)^{w{z)). 

Following [Kap] (sec also the earlier work [Ble]), we want to find an analytical 
matrix-valued function iy(z) (gauge matrix) and an analytical function w{z) of 
the change of variable such that (7.2) reduces to the following model equation 
for $(u;): 

(7.18) ^'{w)=n(^1 l^^w). 

This model equation is nothing else than the Airy equation ip"{w) = N'^w'4){w) 
written in the matrix form. The functions W{z) and w{z) are determined 
itcratively in powers of and we consider approximation of the order of 
A^~^. In this approximation the functions w{z) and W{z) are defined as follows. 
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1. Function of change of variable. Put 
(7.19) w{z;z2)= l^J^^i^{z)dz] 



where 

(7.20) i^{z) 



,0 9^ 
2 



1/2 



and Z2 is a root of ^{z) close to 2:2, 

K^2^) = 0, z^-Z2 = 0{N-'). 

The argument Z2 in w{z; Z2) reminds us that this function is constructed in a 
neighborhood of the turning point ^2- Observe that the function i''^{z) can be 
written in terms of the elements of the matrix A^(z) as 

(7.200 Az) = -detA\z) + N-' f (a?i)'(^) - a?i(^)^#^) 

which coincides (up to terms of the order of N~'^) with the potential of the 
Schrodinger equation equivalent to (7.2) (cf. (3.26,7)). In (7.19) and (7.20) 

we take the branches of the power functions that are positive on the positive 
half-axis. Since z^ — Z2 = 0{N~^), for large N the function w{z; Z2) is analytic 
for z & 0,2 and 

w'{z; Z2) / 0, z & 0.2- 



2. Gauge matrix. Define 



1 



(7-21) W{z-Z2) = ^:^\_alM ^ 

where 



w'{z; Z2) 



> 0, z> Z2 



N 



The function W{z; Z2) is analytic in z e ^l2 for large A'". 

To get a desired multiplicative jump matrix S between Of and ^2 we use 
different solutions to the model equation (7.18) in the domains fif ^^'^ ^2- 
They are defined as follows. 
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3. Model canonical solutions. Define 
(7.22) 



where 

(7.23) y^{z)=ki{z), 

Here Ai (z) is tlie Airy function, i.e. the solution of the model equation 
(7.24) 

1 / 22'^/2 \ 
Ai"{z) = zAi{z); Ai{z) ^ ^^^xj^ ^^P I §— 1 ' ^ ^ go. 

It satisfies the relation 

Ai (z) + e^^^/^ Ai +e-2'^*/3Ai (e-'^'^^l^^ =0; 

hence 

y\{z) - y2{z) = -iyo{z), 
and this implies the equation 

(7.25) $„(^) = $,(^)(^^ 

Now we define the turning point functions 

(7.26) {z; Z2) = CW{z; Z2)$m(^(^; ^2); ^2), 

where C is a constant which will be chosen later. 
In Hi we define 

(7.27) {z; zi) = CW W(^; 2i)$„,d(^«(z; Zi); Zi), 

where C^^^ is a constant to be chosen later, and w{z;zi), W{z;zi), and 
^u,d{z;zi) are defined as follows. The function of change of variable 

^J\{z)dzj , 
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where ^{zi ) = 0, zi — z^ = 0{N ^). Observe that v{z) < when < z < Zi , 
hence ^ 

/ i'{z)dz>0, 0<z<Zi, 

and we take the branch of ^^/^ in (7.27') which is positive for ^ > 0. The 
model solutions ^u,d{z',zi) are defined as 

(7.27") 

<^u{z,z,)-\^ N-^' ) [y'o{N'/'z) -y',{N'/' z) ) [ {2n) 



1/2 , 



<id{z;zi) 



N'"' \(y,{N^'^z) -y,{N^'^z)\({2'Kr^'^ 
iV-i/« M j/i(iV2/3^) -j/i(iV2/3^) M (27r)i/2 



Comparing this with the model solutions (7.22) in ^2, one can notice that we 
change y\ for —y2 in the formula for and 2/2 for —yi in the formula for 
Note that equation (7.25) still holds. The gauge matrix W{z; zi) is defined as 



(7.27'") W{z-zi] 



w'{z;zi) \ — 1 



1 






w'(z;zi) j 




a'l,(z) / 




/ aUz) 


/ w'{z;zi) 



■ > 0, < z < , 
and it is analytic in z e fli for large N. Now we define 

(7.28) *0(,) = P^P^-;-i,2), if ^€f^I'2^ 

1 (-i)V3*0(-^)a3, if ze{-n). 

We reiterate that our use of equation (7.2) has been purely formal: we 
only refer to it as a motivation for the introduction of the function '^'^{z) by 
explicit formulae (7.9) and (7.28). 

Properties of"i!^\z). As defined by (7.9) and (7.28), the function '^^{z) is 
a piecewise analytic function in C and has the following properties. By (7.11) 
and (7.28), 

(7.29) *°(z) = (-I)"a3*°(-z)a3. 
By (7.25), 

(7.30) *°(z) = *°(z)5, zenn{lmz = 0}; S = (j^ "l""') " 

Now we want to evaluate the jump matrices on di^i and dfl2- Observe that the 
WKB function (7.9) has two free constants Cq and Ci, and each of the turning 
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point solutions (7.26), (7.27) has one free constant, C and C*^^-*, respectively. 
We show that we can adjust these constants to make the jump matrices close 
to I. 

Lemma 7.2. // 

(7.31) C7o = -^, C, = -'^- CW = (-1)^C, k=[l 

then 

(7.32) *° (z) = (/ + 0{N-^)) *o (z), 
uniformly with respect to 

z G U 8^2 U {-d^i) U {-8^2) , 

where 

-dQj = {z: -zed^j}, j = 1,2. 



The proof of this lemma is given in Appendix C below. Consider now the 
asymptotics of "^^{z) at infinity. According to (7.15) and (7.31), 



(7.33) ^°{z) 
with 



/y. e-(^-"ln^+iV7„+l^)a3 



^""Vo ij' ^'~[iKy/^ 

We want to adjust the asymptotics to the one (5.16), (5.17) in the Riemann- 
Hilbert problem. Define 

(7.34) A° = iV7n + + C = (i?°)-V4(2;r)V2. 
Then (7.33) is rewritten as 

(7.35) *0(z) ~ ff; e-(^-i-+^^)-3, 



\k=0 



with some coefficients such that 



(7.36) r° - ( J j > r? - f ^^0^^1/2 J 1 . 
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Proof of Theorem 1.1. We want to prove that "^^(z), defined in (7.9) and 
(7.28) with C, Co, and Ci given by (7.31), (7.34), is close to ^'„(2). Observe 
that by (5.16) and (7.35), the hmit 

(7.37) r = ^Inn ^i>n{z) [*°(^)] = roe^^"-^")^^ (rg)-^ 
exists and detF / 0. Denote 

(7.38) X(z) = r-i^„(z)[*0(z)]-'. 
Then 

(7.39) hm X{z) = I. 

2— »0O 

The multiplicative jump matrices of "^niz) and "^^{z) on $1 n {Imz = 0} 
coincide (both are equal to S); hence X{z) has no jump on H {Imz = 0}. 
On j = 1, 2, the function ^niz) has no jump; hence 

x+{z) = r-i*„(z) [^liz)] = r-'^n{z) [{I + o(iv-i)) (^)] -1 

= T~'^n{z) (2)] (/ + 0{N-'))-' = X^iz) {I + 0{N-'))-' , 

so that 

(7.40) X+{z) = X_{z) {I + 0{N-'^)) , if zedniUdn2. 
Let 

(7.41) L = {z: lmz = 0, z ^ Q U (-Q) } . 
On L, the function *°(z) has no jump; hence 

(7.42) 

x+{z) = r-'^n+{z) [*°(2)]"' = r-^^n-{z)s [*°(^)]"' 

= r-'^n-{z) [*°(-2)]~' ^''{z)S [*°(z)]"' = X_{z)Q{z), 

where 

Q{z) = ^\z)S [^\z)]~\ zeL. 
When z e L, *°(z) = *wkb(-2;); hence by (7.9), 

Q{z) = ro(z)e[-^«(^)+^il'^3Se[^«W-^il^3To-H^) 

^l -27rie-2^«W+2Ci \ 
= W(i 1 )To-\z) 

/n _ie-2JV€(2)\ 
= I + Mz)[l )T,-\z). 

Since detro(z) = 1, we obtain from (7.10) that \To{z)\, \Tq^{z)\ < C; hence 

g(z) = / + o(e-2^^^«^)) 
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and 

(7.43) X+{z) = X_{z) + O (^e-^^'^^^^^))) , z e L. 

Thus, we have the fohowing theorem. 

Theorem 7.3. The function X{z) = r-i*„(2;) satisfies 

hm X(z) = I 

2— >0O 

and 
(7.44) 

X+{z) = X^{z) {I + 0{N-^)) , if z e dniUdn2l^{-dni)U{-dn2), 

X+(z) = X_(z)(l + o(e-2^I^e«W)), if zeL; 

i. e. , there exists a constant C > independent of N and z such that 

[x^{z)]~^ x+{z) - 1 <CN-^, if z G u u (-50i) u (-5O2), 

[X_{z)]-^ X+{z) - I <Ce-^^^^^^'\ if zeL. 

Remark. Observe that there exists cq > such that Re^(z) > 
co(l + \z\)^ > for all 2 G L (cf. (7.7')); therefore the error term on L is 
exponentially small as N ^ 00 and also as z ^ ztoo. 

By standard methods (see e.g. [DZ2], [DIZ] and also Appendix D below), 
Theorem 7.3 yields the asymptotic equation, 

(7.45) xW = / + o(^^^^), .ec, 



(7.46) ^Jz)=r{I + 0{ .„ , „ I \<S"(z). 



so that 

= r ( T -i-n I 

^N{l + \z 
This implies that 

Equating the coefficients at and at z~^, we obtain 

(7.47) ro = rrOe^^3^ ^^^^_^o^ 

Ti = r [r? + o(A^-^)rO] e^^K 

Excluding T we obtain 

e'^'''{To)-^Tie-^^^ = (r[])-^r? + 0(iV-^) 
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We have 



hence 



Thus, 





Rnc-^^ 



e2'5' 
i?„e-2« 

which imphes that 

(7.47') Rn = Rl + 0{N-'), A„ = A°+0(iV-i). 

In view of the definition of (see (7.1)) the first equation yields (1.14). By 

(7.47) , 

r = r„.-»'.(rgr' = (^; ^„,^,_°,(^,„,) =/ + o(iv-'); 

hence by (7.46), 

(7.48) *„(z) = r (/ + o ) r = / + O (7V-1) . 

All the statements of our main Theorem 1.1 follow directly from equation 
(7.48). In what follows we are presenting the corresponding detailed deriva- 
tions. 

To derive the scmiclassical asymptotics (1.15) of the function tpniz) be- 
tween the turning points zi and Z2 we use the following lemma. Let 5 > be 
an arbitrary fixed number. Define 

Qq = Qn{z: zi + S <Re z < Z2- 6} 

and 

n^ = non{z: lm^>0}, n^ = nor]{z: Im2;<0}. 
Lemma 7.4. The function '^^{z) has the following asymptotics in ^q: 

"12 (2) "12(2) ^ 



where ^{z) is defined as in (7.6), (7.7'), 

(7.49') ^(.) = -^n( ''W 

2 V 0^2 (^) 
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and 



= Qi' 



i 1 
-i 



/(27r)-V2 
(27r)V2 



Proof of this lemma is given in Appendix C below. Due to (7.48), we can 
replace with *n(^) in (7.49), which gives the semiclassical asymptotics 

of *n(2) in ^0- Thus, we have the following theorem. 



Theorem 7.5. In ^o, 

(7.50) *„(.). (^+^(^'^))^?2(^)^^^^ 



(i20)V4 V2m(^ 



1 







^ g(iV^(.)+<i(.)).3Q(2)p^ i/ 

Remark. Let ^0(^)1 Co (-2), c?o(-2) denote the analytic continuation of the 
functions lJ.{z), S,{z), d{z) from to across the interval [zi + d,Z2 — S]. 
(Note that fio{z) = —/^{z), S,oiz) = —^{z), and do{z) = —d{z) for z G Oq-) 
Taking into account that by (3.14) the function il)n {z) is the element 11 of the 
matrix- valued function '^^{z), we derive from (7.50) that 



(7.50') Vn(^)= (/ + 0(iV-i)) 



a\2{z) 
(27r)V2(i?0)i/4 



1/2 



je^«o(2)+<io(2) + ^-{N(,o<,z)+do{z)) 



if z e fto. 



The function N^q{z) + (io(^;) is pure imaginary for real z in the interval 
zi + 5 < z < Z2 — 5 (cf. (7.7') and (7.49')); hence, when we restrict the 
asymptotics (7.50') to the interval zi + S < z < Z2 — S, the equation. 



(27r)i/2(i?0)i/4 ^2/i(z 



le 



+ e 



-{Nio{z)+do{z)) 



+ o(iv-i) 



follows. Since 

we obtain that for zi + (5 < z < 2:2 — (5, 



(7.51) 



'ipn{z) 



0FVAy (l-x^ 



!)l/4 



cos (<^{z) + 0{N-^] 
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where ^{z) = i ^ [N£,o{z) + do{z)]. To finish the proof of (1.15) we use the 
following lemma. 



(7.52) 



Lemma 7.6. If zi + 6 < z < Z2 — 6 then 

(n + i) / sin 20 



where (p = arccos q and x = arccos r and 



gz' + t 



2y/Yg-tq 



X' 



(-i)"x 



(n + i) 



Proof of this lemma is given in Appendix C below. Since 



we obtain 



Vl-x^ = sine/) (1 + 0{N-^)) ; 
hence for zi + 5 < z < Z2 — S, 



Ipniz) 



1 /5\i/4 ^ 



x/tF \\J ^sin( 



cos 



sin 2(j) 



4 4 



+ O(Ar-i) 



Formula (1-15) is proved. 

To prove formulae (1.17), (1.18) we use the following theorem. 

Theorem 7.7. The function '^n{z) has the following asymptotics: in 
C\ (OU (-0)), 



(7.53) ^n{z) 



«12 \ 



'-12 



\ 0- ^ 

where the matrix Cn does not depend on z, and 

c„ = / + o(i), c:w=/ + o(^j^). 

In Qj, j = 1,2, 

(7.54) ^n{z) = + J)^^(-ir ^.); ^ . ) , 



(i?0)l/4 



^Gi^J'^ i = l,2, 
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where ao = (2 - j) [§] , 

(7.55) w{z;zj) = ^ J ^u{z)dzj 



1/2 



and the gauge matrix W{z;zj) and the model canonical solutions ^u,diz',Zj) 
are as defined in (7.21), (7.27'") and (7.22), (7.27") respectively. 

The proof of the theorem fohows from the definition of '^^{z) and (7.48). 

Theorems 7.5 and 7.7 provide semiclassical asymptotics of ^n(^) in the 
whole complex plane. The formulae (1.17) are derived from (7.53) in the same 
way we derived (1.15) from Theorem 7.5. Let us show that (1.18) follows from 
(7.54). For the sake of definiteness consider j = 2. Restricting (7.54) to the 
element 11, we obtain that if z G then 



(27r)V2 
(i?0)V4 \w'{z;z2 



al^(z) 



1/2 



(1 + £i)7VV6 Ai (w;(z; Z2))(27r)-V2 
+ e2N-^/^Ai'{w{z;z2))\ , 



where £1,2 = 0{N ^). This gives 

N^^^z 



[(1 + n) Ai {w{z- Z2)) + r-2 Ai' {w{z- Z2))] , 



{w'{z;z2)y/^ 

where ri = 0{N~^) and r2 = 0{N^'^/^). Formula (1.18) is proved. 
Finally, by (4.14), (7.47'), (7.34), and (7.14), 

hn = exp (2An) = exp (2A° + 0{N-^)) 

= exp ( 2iV7„ + ln27r + + 0{N-^)\ 



Nt^ NX 



9 



1 + In ^ ) + 0{N-') 



L 45 2 

which proves (1.23). Theorem 1.1 is proved. 

Appendix A. Proof of formula (7.8) 

By (7.6), 
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and 



diag T-^T' = 
In our case = —ai2'^ hence 



V 







"12 I "21 



(A.1) 
and 

Now, 



/.2 = (a50'-(a?2)' 



diag T-^T' = - 



7O \'„o /I 
h2 \ "12 ' 



"12 \ "12 



H-a1^J 2fi \0 1 



'12 



By (A.l), 



«?2(«?2)' = «ii(«ii)'-w'; 



hence 



«?2 A' «?2 _ [<{<)' - - - - f^'] if^' - 



'^0 /-^O \/,,2 /„0 \2,,/ 



(/"-ail) M 



X ' 



= In 



A* -Oil 



Thus, 



and 



hence 



diagT" 



-V = lfln^i^V/ 



r(z)=/ dl^gT-\u)T'{u) du = \\n ^^^^ 



U=Z2 



-1r{z) ^ ^ /x(^) -aii(^) j 



Formula (7.8) is proved. 
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Appendix B. Proof of Lemma 7.1 



Note that 



(B.l) 

— all ^ + Vx"^ — i + c 



2fi 



2Vx^ - 1 



xo = - 



2V¥ 



> 1. 



Let ri* denote the region, 

n* = {z eC: Rez > 0} \ [2;i,^;2], 
and define the analytic functions on Q* by the equations, 



(B.2) 

and 
(B.3) 



Ci{z) = 2Vx2 - 1, C2iz) =x + Vx'^ - I; 



(^{z) = x + Vx^ — 1 + c . 



X 



t + gz' 



These functions generate the conformal mappings, z ^ ^, of J7* onto the 
regions of the complex plane ( depicted in Figure 3 below. 

Let us also denote by C^/^ the branch of the square root of C defined on 
the complex plane C cut along the ray (— oo, 0] and fixed by the inequality. 



Note that 

(B.4) 

and 

(B.5) 



— TT < argC < TT. 



\ 1/2 



2/i 



From Figure 3 it follows that 

(i) The functions (-^ ^^^{z), (2 ^{z)i Cs^'^i^) single- valued on 

n*\[o,zi]. 

(ii) The jumps of the functions (^^ ^^'^(z), (2 ^(-2)) and Cs^'^i^) across the 
interval [0, Zi] are described by the equations. 



(B.6) 



n n 

cl^\z) , = i-ir+\i/\z) 
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Fiszure 3. The conformal mappinas related to ^ (z) 

WKB 



(iii) The product, C3 (•2)Ci (-2), being analytically extended to 
C \ [—Z2, Z2] across the imaginary line becomes an even function. 

n 

(iv) The function, C2 ^(-z), being analytically extended to C \ [—-22) -^2] 
across the imaginary line satisfies the symmetry equation, 

n n_ 

(B.7) ^p^-,) = ^-l)n^-^,). 

Prom the properties (i)-(iv). Lemma 7.1 follows. 
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Appendix C. Proofs of Lemmas 7.2, 7.4, and 7.6 

Proof of Lemma 7.2. It is easy to see that the function defined by 

(7.9) and (7.28) satisfies the equation, 



(C.l) ^"(z) = ^o(z). 

Taking also into account the z —z symmetry (cf. (7.29)), we conclude that 
it is enough to prove the upper-right part of relation (7.32). We shall start 
with discussing the jump over the arc (see Fig. 2). 

For every e > 0, in the sector — tt + £ < arg z < tt — e, the Airy function 
is written as 

(^) = ^ / ,/ X exp[-?7(2;)] , -TT + e < arg z < TT - £, 
2A/7r?7'(z) 

where 77(2;) is an analytic function, with the asymptotics 
(C.2) 



r?(2) = ^z3/2(l + 0(|z|-3)), ,7'(2) = zV2(i + o(|^|-3))^ ^^^^ 

(cf [Ble]). Hence 
(C.3) 

"'"^ i^'^'') - (-^) + ^ (^")) ' 

iV-V6 Ai' (iv2/3,) = -|^exp (^-^^^ (1 + (iV-^)) , 

|z| > £, — TT + £ < argZ < TT — £, iV — GO. 

By (7.23) this implies that 
(0.4) 



z| > £, — TT + £ < argz < TT — £, N ^ oo, 
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and 
(C.5) 



27r 



|z| > e, — TT + £ < arg z — <7r — e, iV— ^ oo. 



Hence, by (7.22), 
(C.6) 



1/4 u/4 (/ + 0(A^ '))exp P, 



2Vvr 



TT 



|z| > e, — — + £ < axgz < TT — £, iV — oo. 



where 

Thus, 
(C.7) 



^ - ' (27r)V2 



1 /t/;-V4(^;^2) 

^u{w{z; Z2); Z2) = [ 1/4. x 1/4. 

2^/K \-w'-'^{z;z2) w^i^{z;z2) 



under the assumption that 



(C.7') 



\w{z;z2)\ > e, 



TT 



— + e < argw{z; Z2) < TT — e, N 



00. 



Observe that this assumption on w{z;z2) holds for all z G C2 when N is 
sufficiently large and 6 > (the vertical half-axis of the ellipse Q) is sufficiently 
small. Indeed, from the definition of the function w{z;z2) given in Section 7 
(see (7.19), (7.20)) it follows that ±w{z;z2) > when ±{z - z^) > 0. In 
addition, w'{z\Z2) 7^ for z G ^2 and w'{z-,Z2) > for real z e ^2- This 
insures that there exist 6 > and e > 0, such that the assumption (C.7') on 
w{z; Z2) holds when z G C2 , for all sufficiently large N. 



238 



PAVEL BLEHER AND ALEXANDER ITS 



By (7.19), 



(C.8) 



■ w 



3/2 



{z;z2)= j ^J ii^{u) + N-^U^^){u) du 



+ 

70 



H'^{u) + N-W(^){u)du, 



where 



and 7o is a contour which starts at the point zq = ^^^^ on the lower side of 
the cut [zi, Z2], goes around the point and ends at the point Zq on the upper 
side of the cut [zi, 22]- The contour 7(2) is a simple curve which goes from the 
point zo to the point z and does not intersect any of the cuts. On 70 U 7(5;), 

(C.9) ^J,,^{u) + N-W){u) = ix{u) + + O (7V-2) , ue^o^^{z)- 



2Nii{u) 



hence 
(C.IO) 



+ 

70 Jl{z) 



+ 



= f\{u)du + N-^ [ 
= + N-^d{z) + O (iV-2) 



fj,{u) du 

1 U^^\u)du 
J Mu) 

U^^\u) du 



+ o (Ar-2) 

+ O (Ar-2) 



where 



du = 2^ — ln(a; + a/x^ — 1) 



2^ ' 



(cf. (7.6), (7.7')) and 
(C.ll) 



d{z) = f 



U'^^\u) du 



1 r Ki( 

2 Jz2 a?2 



«12(^) 



+ 



(a?, (n))'a02(^)-a;i(n) (a°2(^))^ 
(n)VKiH)^-K2WP 

«12(«) 



- 1 



1/2- 



= "2 Ina(^), 
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where 



a{z) = 
Note also that 
(C.12) a{z) = 

where 



1 /A 



z Vf/ Vc 



Substituting (C.IO) into (C.7), we obtain the following asymptotic relation: 
(C.13) 

1 (w-y\z-Z2) w-y\z-Z2)\ 

^u{w{z-Z,)-Z,) = ^^_^,/,^^_^^^ ..V4(,.,,) J 

X (7 + O (iV-i)) e-(^«(-)+<^W)^3p^ 

(assuming again (C.7')). Let us consider the function (z; Z2) for 2; in a small 
neighborhood of the arc C2 • In this case, as we discussed above, the condition 
(C.7') holds (for all sufficiently large A''); hence we can use (C.13). We want 
to check that for z e C2 the turning point solution *?pp {z; Z2) coincides, up to 
terms of the order of N"^, with the WKB solution (7.9), with some constants 
Cq,C\. Differentiation of (C.IO) gives that 

(C.14) w'{z- Z2)w^''^{z- Z2) = ^'{z) + O(Ar-i) ^ ^ O(iV-i). 

Hence we obtain from (7.21) that 

(C.15) 

'aO,(.)xV2 



Wiz;z2) 



H{z) 




and by (C.13), 



W{z-Z2)^u{w{z-Z2)\Z2) = :—j= 



1 (aUzWi'i ,! ^ ^ U 1 1 




X (/ + 0(Ar-i))e-(^?W+<^W)^3p 

20F V /i(^) / V«"H^) a(^) 
X (I + 0(iV-^))e-(^«(^)+'^(^))'"3p_ 



-1 1 
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To adjust this formula to the WKB solution (7.9), we rewrite it as 



In virtue of (C.ll), this gives that 

(C.16) 

W{z;z2)^u{w{z;z2);z2) 

1 U{z) - aUz) \^' ( 1 --'^^Aii + oiN-')) 
"V^l ) 1 )^I + 0{N )) 

V 27r 



(of. (7.9)). Thus, by (7.26), 

(C.17) *Tp(^;^2) = ^ro(z)(7 + 0(Ar-i))e-^^(^)^3-(i/2)(in2^)a3^ 



uniformly for z in a small neighborhood of the arc d^- Therefore, if we take 
in (7.9) 

(C.18) Co = -^, Ci = -(l/2)(ln27r), 



then 

(C.19) v^J^p {z- Z2) [^wkb(^)]~' = / + O(iV-i) , z G Cf. 

Because of symmetry relations (C.l), we obtain simultaneously that 
(C.20) {z- Z2) [^'wkb(^)]"' = / + O(iV-i) , z e Ci 

Consider now the jump across the arc C". By (7.23), 



(C.21) y2{z) = y^{z)- 

hence by (C.5), 

(C.22) 

ivV6y,(iv2/3,) = _A_ (^-^^ (1 + O(iV-i)), 

N-y%'2{N'l'z) = l^exp (^^^) (l + 0(iV-i)), 

27r 

1^1 > £, — TT + £ < arg2; + — < TT — £, A' ^ oo. 

o 
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Combining this formula with (C.5) and recahing the definition of the function 
^u{z;zi) in (7.27) we have 

(C.23) 

II 7!" 

|2;| > e, — TT + e < argz < — — £, N ^ oo, 

where P is the same as in (C.6). Thus (cf. (C.7)), 
(C.24) 

^ / / XX ^ f w~^/'^(z: Zl) —w''^/^(z:zi)\ ,^ ^ /,r-i\\ 

*„K.;.,);^i) = 2^ _„i/4(L/,,) j(^ + o(iv ■)) 

provided that 

vr 

(C.24') \w[z\z\)\> — TT + £ < argu)(2;; Zl) < — — e, — > oo. 

From the definition of the function w{z\z\) given in Section 7 (see (7.27')) it 
follows that ibti;(z;zi) > when ib(z — z^) < 0. In addition, w\z\z\) ^ 
for z G rii and i(;'(2:;zi) < for real z G Oi. This secures that there exist 
6 > and £ > such that the assumption (C.24') holds when z G C", for all 
sufficiently large N . 
By (7.27'), 

(C.25) -u;3/2(z;zi) = T \U\u) + N-^lJ^^){u) du 

= / y^/i2(u) +iV-iC/(i)(^x)dn 

\ I J ix'^{u) + N-W){u) du, 
2 Jt-, 



+ 

•} I 

71 



where 71 is a loop which goes clockwise around the interval [zi , ] • 
loop 71, the integrand of the last integral can be estimated as in (C.9). Hence 

(C.26) 

\ I Jfi^{u)+N-W){u)du= r\{u)du + N-' rE^p^ 

^ J-yi J Zl J Zl lixyu) 

+ 0(7V-2) 

= -e(^i) - iV-id(zi) + o (iv-2) , 
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where the functions ^(z) and d{z) are the same as in (C.IO), (C.12) and their 
values at zi are taken from the upper half-plane. Using the notation introduced 
in Appendix B, we can rewrite the functions ^{z) and d{z) in the form 



n 

2N 



x-s/ — 1 — In C,2{z) and d{z 



This yields (cf. Fig. 3) that 



Z7rn 



and d{zi) 



2 IC3(2l)| 



-llnil, 

2 z 



where = 1 if n is even, and = if n is odd. The last equations together 
with (C.25) and (C.26) lead to the following representation of the function 
w{z; z\) in the neighborhood of the arc (cf. (C.IO)), 



(C.27) 



t^^/\z,z,)=i{z) + N-^d{z) + ^ + '-^ 



+ 0{N-^). 



Substituting (C.27) into (C.24) and assuming (C.24'), we arrive at the following 
asymptotic relation (cf. (C.13)): 

(C.28) 

(-1)^ ( w-'^l^{z;zi) -w-'^l\z;zx)' 
■w^l^{z\z\) —w^l'^{z\z\) 



^u{w{z\zx)\zx] 



20F 



(/ + o(iv-i)) 



k = 



As discussed above, for z in a small neighborhood of the arc Cf and suf- 
ficiently large AT, condition (C.24') holds; hence we can use the asymptotic 
formula (C.28) in equation (7.27) for (-^S -^i)- compare 'I'^^P (-^' -^i) with 
^wkb(-2)) we also need (as in the previous case of the arc C^) the asymptotics 
of the gauge factor W{z\ z\) in (7.27). Similar to the previous case (cf. (C.14)), 
differentiation of (C.27) gives that 

vJ{z- zr)w'l'^{z- zi) = ^{z) + O(Ar-i) ^ ^(^^ ^ O(Ar-i), 
which in turn implies the relation, 

(C.29) {w\z- z{)fl^w^'\z- zi) = -tx^'\z) + 0{N-^). 

(The appearance of the sign "— " on the left hand is due to the condition 
i(u;'(z;zi))V2 > 0, < z < ; see (7.27'").) Prom (C.29) and (7.27"') it 
follows (cf. (C.15)) that 

(C.30) 

W{z-zi)=i 




w^/'^{z] Zl) 




w 



-V4( 



Z\Zl 



J(/ + 0(A-i)), 
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and by (C.28), 

W{z;zi)^u{w{z;zi);zi) 

i-li nil- - 

X (/ + C>(iV-i))e-^^3e-{ive{^)+dW)a3p_ 

ITT 

Since e""*^^ = — iua, we obtain 

(-I)'' /a° {z)\^^'^ ( ^ ^ \ / 1 -1 



1 




\^ (/ + 0(Ar-i))e-(^«(^)+<^(^))'^3p 



X (/ + C)(Ar-i))e-(^«^)+<^))^3p 



2^ V Ai(^) y Va~n^) a(z) 

The rest of the derivation is exactly the same as in the case of the arc C2 , and 
it yields the following formula (cf. (C.16)): 

^^ i i ^ ^ (-1) VM^) " ^ '^V 1 



/27r 

Thus, by (7.27), 



(C.31) *|p(z;zi) = ^^^^^IJ^ ro(z)(/ + 0(Ar-i))e-^«(^)^3-(i/2)(in27r)a3^ 

v 27r 



uniformly for z in a small neighborhood of the arc Cf (cf. (C.17)). 

To ensure the same WKB-asymptotics both for *?pp {z; Z2) and (z; zi) 
we take C^^) = (-l)^^. Then 

(C.32) ^-^'^ (z; zi) [*wkb(^)]-' = / + O(iV-i), z G C^^'- 
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Our last step is to prove that the turning point solutions ^'^p {z; zi) and 
^'^p {z; Z2) coincide, up to the terms of the order of , on Iq = dVL\r\d^2 (see 
Fig. 2). By (C.19), (C.20) and (C.32), -^V^-p {z; Z]) and -^V^-p {z; Z2) coincide, up 
to the terms of the order of A^~^, at the end point of Zq, z = z^ + ih. Since 6 > 
can be chosen as small as we want, this is true for all 2: = + ic, < c < h. 
The problem is to get a uniform estimate of the error term up to c = 0. We 
will consider a more general problem. 

Let (5 > be an arbitrary fixed number such that 5 < {z2 — Zi)l2 and let 

(C.33) ^o = ^r\{zi+5 <Kez < Z2- 6}. 

Then Zq C Oq. We will prove that 

(C.34) *^p(z;zi)[^|p(z;z2)]-' =/ + 0(iV-i), z G J^o, 

and 

^=I + 0{N-^), zG^lo- 

Observe that w{z; Z2) is analytic in z G (On {Re z> zi+ 5}) if N is sufficiently 
large, and w'{z;z2) 7^ in this region. Hence the ^'^p(z;z2) are analytic 
there as well. Similarly, the '^^^ (z; zi) are analytic in the region n {Rez < 
Z2 — 5}. Thus, all the functions ^'^p (z; 21^2) are analytic in Qq- Because of the 
symmetry equation (C.l), it is enough to prove (C.34) only. The idea of the 
proof of (C.34) is to find a WKB function ^'wkb('^) ^^'^^ that for A; = 1, 2, 

(C.36) ^>^^p{z■Zk)W^K^{z)\-^ = I + 0{N-^), zeno. 

To construct the function ^wkb(^) consider an auxiliary turning 

point solution 

(C.37) {z; Z2) = CW{z; Z2)Mw{z- Z2)), 

where 
(C 38) 

yi(iV2/3^) y2(A^2/'^) W (27r)-V2 Q \ 

y[{N^/^z) yi,{Ny^z) )[ (27r)V2 ) ■ 

The key point here is that in contrast to ^u,d{z',Z2) in (7.22), the matrix 
elements of ^o{z) are defined in terms of yi{z), y2{z) only (and not yo{z); cf. 
(7.23)), and this allows us to find the asymptotics of "^^p {z; Z2) as N ^ 00 
for all z & Uq. The function y2{z) is an entire function. Changing the branch 



(C.35) 



'^ipiz;zi) *Tp(-z;22 



Mz) 





7V-V6 
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of the root functions on the right-hand side of (C.22) we obtain 
(C.39) 

II 47r ,^ 

\z\ > £, — TT + £ < arg z — < TT — e, iV— ^ oo. 

Combining this formula with (C.5) we obtain 

, , TT 57r 

|2;| > e, — + £ < argz < — e, N ^ oo. 

Let ^oiz), l^oiz), and dQ{z) be the analytic continuation of the func- 
tions ^{z), n{z), and d{z), respectively, from Oq to Qq, across the interval 
[zi + S,Z2 — S] (see equations (7.4), (7.6), and (C.ll), (C.12)). Then equations 
(C.IO), (C.14), and (C.15), with ^{z), iJ,{z), and d{z) replaced by S,oiz), Hq{z), 
and do (2;), hold for all z € ilo- In addition, there exist 6 > (the vertical 
half-axis of O) and £ > such that for all z in Q,q and large N the inequalities 

(C.41) \w{z]Z2)\> e, < axgw{z;z2) <^ - e, 

hold and hence one can use the asymptotics (C.40) for the function ^q{w{z;z2))- 
All this together leads to the asymptotic equation, 

(C.42) 

X (/ + 0{N-^)) e(^io{z)+do{z))a3p^ ^ g 

(cf. (C.16)). Since iyiiz) — iy2{z) = yo{z), we obtain that 

(C.43) ^.iz;z2) = Mz)P-'(\ 1)P, ^=(^'T'' (2.V)' 
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hence 
(C.44) 

{z; Z2) = CW{z; Z2)^uiwiz; Z2)] Z2) 



Thus, if we define 
(C.45) 



then 

(C.46) {z; Z2) = {I + 0{N-^)) *Wb(^)- 

Consider now in fio the turning point solution '^'yy{z\Z\). For z in J7o 
the same equations (C.27), (C.29), and (C.30), with the functions ^(2), 
d{z) replaced by the functions £,o{z), iio{z), dQ{z), hold. Also, there exist 6 > 
and £ > such that for all z G Qq and sufficiently large w(z; zi) lies in the 
sector 



(C.47) -^+s<&igw{z-zi) <~-e 



and |u)(2;;2;i)| > e. This implies that 



W{z;z,)Mw{z;z,)) = -^{-l)''(^Y' ( Xi^) Mol)V/ I 

X (/ + 0(iV~l))e~^'^3g-(Af?o(2;)+rfo(2:))o-3 p 
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Since e 2 '^s _ we obtain that 

(C.48) 

W{z;zi)^o{w{z;zi)) 



»12(^) »12{^) 



X 



_|_ (9(jV-l))g-(^5o(z)+(^o(z))o-3 p 



(-1)' /a?2(2)V'V ' " Vl i 



Taking into account that (cf. (C.43)) 

we derive from (C.48) the asymptotic equation 
^-Ip {z;zi) = {-l fCW{z;zi)Mw{z;zi)) 

C (aUz)V^'( ^ M^l 



-l\ f i 
-1 j I -i -1 



1/2/1 



In other words, 

(^; zi) = (7 + O(Ar-i)) *^KB 

hence 

^-Ip (z; zi) [*|p (z; Z2)]"' = / + 0(iV-^), z G J^o- 
Lemma 7.2 is proved. 
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Proof of Lemma 7.4. This lemma follows directly from the formulae (C.45), 
(C.46) (and symmetry (C.l)) if we substitute into these formulae the value of 
C given in (7.34). 

Proof of Lemma 7.6. By (C.8), (C.IO), 
(C.49) 

^o{z)+N-'do{z) = \l [ +[ I ^J^^u) + iV-iC/(i) (u) du + (7V-2) , 
[ ^ -'70 -J-yiz) ) 



where 



4 N ' 



/x^(z) = 



C/(^)(.) = (a?,(.))' - a?,(.)i%f = + I + 9Rl 



so that 
(C.50) 



li^{z) + N-^U^^\z) = 



{gz^ + tf [n + i) g 



N 



= z 



-X'g 



N- 



{-ir^t^-A\g 



= ni{z) - N~'a^ 



where 



fii{z) = X'gz\q''-l), q = 



gz' + t 

2v^ 



{-lY^f^-AXg 



Hence, 



[ ^ -'70 -'7(z) J 



= r (Mu)-N-'-^) du + 0{N-') , 
where z^q is a close-to- Z2-zero of Hi{z). Now, 

/ ijLi{u) du = \/X/g / — 1 ud-u = is/yTg \ \/\ — q 

] j .N I 

^ Znn Ze^n Znn 



udu. 
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The substitution v = reduces this integral to 

iX' 



/ du = iX' / \/l — v^dv = — yqy/l — q^ — arccos qj 

Jz^Q Ji 

iX' ( sin 20 



2 

Similar computation gives 



arccos q. 



du . 2du 



Mu) Jz^^ Uy/iX'g - {gv? + 



I arccos r ^ ix 



2y/Yg-tq 

r = — ; , Y = arccos r. 

2^fX!gq-t ^ 



^ [N^o{z) +doiz)] - ' 



Thus, 



2 

and Lemma 7.6 is proved. 

Appendix D. Proof of estimate (7.45) 

Denote 

s = L u u 5J^2 u {-dni) u {-dn2). 

Theorem 7.3 implies that the function X{z) solves the Riemann-Hilbert prob- 
lem on the contour S, 

(D.l) X{oo) = lim X{z) = I, 

(D.2) X+{z) = X^{z)G{z), zeE, 

with the jump matrix G{z) given by the equations, 
(D.3) 

= *° (z) (z)] ^ if zedniUdn2Lli-dni)U{-dn2), 

G{z) = ~^ if z e L, 

and satisfying the uniform estimates, 
(D.4) 

\i -G~^{z)\<CN-^, if zedVLi\jdn2yj{-dni)u{-dn2), 

\I -G-^{z)\<Ce-'^^^^^^'\ if z&L. 
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The Riemann-Hilbert problem is depicted in Figure 4. 



Figure 4. The Riemann-Hilbert problem for the function X(z) and the contour E 

Similar to the asymptotic problems considered in [DZ2] and [DIZ], the 
crucial fact, which follows from (D.4) and which is indicated in Figure 4, is 
that 

(D.5) \I-G-'iz)\ = 0{N-'), 

uniformly for z G S (compare to the similar analysis in [DZ2] and [DIZ]). 
Simultaneously, inequalities (D.4) yield the corresponding L2-norm estimate: 

(D-6) ||7-G-i|U,(s) = 0(iV-i)- 

Remark D.l. The function "^^{z) is defined by explicit formulae (7.9) and 
(7.28) which, in particular, imply that the functions *±(z) and, in virtue of 
(D.3), the jump matrix G{z) can be analytically continued in a small neigh- 
borhood of dQi U 8^2 U {—dQi) U {—dQ2)- Moreover, in this neighborhood the 
first inequality in (D.4) holds. Indeed, the inequality follows from Lemma 7.2, 
whose proof is based on the estimates (C.17), (C.31), and (C.34), (C.35) taking 
place in the small neighborhoods of the corresponding pieces of the contour 
dUi U • When z e L, 

G{z)=l + T^{z)\^^ '\ jTo-^z), 

(cf. (7.42), (7.43)) and hence G{z) is analytic and satisfies the second inequal- 
ity in (D.4) in a small neighborhood of L, which may include the narrow 
sectors (the angle less than vr/8) about the intervals (—00,-2:2] and [22,00)- 
The indicated properties of the matrix G{z) mean that in the setting of the 
Riemann-Hilbert problem (D.l), (D.2) the choice of contour is not rigid; one 
can always slightly deform the contour S and slightly rotate its infinite parts. 

By standard techniques in Riemann-Hilbert theory (see e.g. [CG]; see also 
[BDT] and [DZ]) the solution X{z) of the Riemann-Hilbert problem (D.l), 
(D.2) is given by the integral representation, 

(D.7) x{z)=I+^.j^pW){l-G-\v))^^, 
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where p{z) = Xj^{z) solves the equation 

(D.8) p{z)=I + C+[p{l-G-^)]{z), zgS, 

and C+ is the corresponding Cauchy operator: 

{C^h){z)= Um f h{ul^_ 

2'-»2,z'6(+)-side7s I/ — r 27rZ 

We note that by (7.38), X{z) is built up from the functions which satisfy 
asymptotic conditions (5.16) and (7.35). Therefore, / — p{z) = O {z~^^ as 
z CO and hence a priori it belongs to the space L2(S). 
Put 

Pq{z) = p{z)-I. 

Then equation (D.8) can be rewritten as an equation in 1/^(2), 
(D.80 {l-B)po = f, 

where the function f{z) and the operator B are defined by the equations, 

/ = C+(/-G-i), 

B:(^^C+[(^{l-G-% <^eL2(S). 

The L2-boundness of the operator C+ (see e.g. [LS]; see also [BDT] and [Zh]) 
and the estimates (D.5) and (D.6) imply that 

I|-B||l2(S)^L2(S) < l|C+llL2(S)^L2(S)l|-f - G^~"^IILoo(S) = 0(iV"^), 

and 

II/IIl2(s) < l|C+|k,(sHL,(E)||/-G-i||i,(s) = 0(iV-^). 

Therefore, equation (D.8') is uniquely solvable in -L2(S) for sufficiently large 
A'', and its solution satisfies the estimate, 

l|po||L2(s) = 0(iV-^), 

or, in terms of p{z), 

\\I - p\\L,i^) = 0{N-^). 
This equation together with (D.5) and (D.7) show that 
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To complete the proof of estimate (7.45) we only need to notice that the 
domain of validity of the asymptotic equation (D.9) can be made the whole 
z-plane because of the indicated flexibility in Remark D.l in the choice of the 
contour S. 

Remark D.2. The kernel in integral equation (D.8) is small (estimate (D.5) 
again), and it is given explicitly in terms of elementary and Airy functions. 
Therefore, equations (D.7)-(D.9) make it possible (cf. [DZ3] and [DIZ]) to 
obtain the full asymptotic expansions for all the three main objects, A„, Rn, 
and Pn{z). 



In this appendix , for the readers convenience, we reproduce with some 
more details the Riemann-Hilbert approach to the orthogonal polynomials sug- 
gested in [FIK]. 

Let si, S3 be complex numbers, |si| + |s3| / 0, and let L denote the cross. 



We assume that the cross L is oriented in a natural way, i.e. from —00 to +00, 
and from —too to +ioo (see Fig. 5). 



Appendix E. The Riemann-Hilbert approach 
to orthogonal polynomials 




iR 






* R 




Figure 5. The Riemann-Hilbert problem for the function Y(z) e NV{z)- 
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Consider the following Riemann-Hilbert problem posed on the cross L for 
2x2 matrix function Y{z) = Yn{z), n G Z+: 

(i) The function Y{z) is analytic on C \L. It has the limits, Y±{z), from 
the lb sides of L which satisfy the jump equation, 

(E.l) 

(1 Se '^^(^) \ NV{z) NV(z) 
Q ^ J=y_(z)e ^"^Se^"^, zeh, 

where 

s = si if 2; G M , s = S3 if z E iR, 

and 

(ii) The function Y{z) is normalized by the asymptotic condition, 
(E.2) lim Y{z)z-''^^ = I. 

z — ^00 

The same arguments as in the case of the Riemann-Hilbert problem (5.16)- 

(5.18) (see Section 5) show that the solution Y{z) of problem (E.l), (E.2) is 
unique (if it exists). In fact, we first notice that since detS* = 1, the scalar 
function detY{z) is a bounded analytic function on C and such that its limit 
as z 00 is 1. Hence, 

(E.3) dety(z) = l. 

Suppose now that Y[z) is another function satisfying (i), (ii). Since the jump 
matrix S is the same for both Y{z) and Y{z), the matrix ratio, Y{z)Y'~^{z), 
is a bounded analytic matrix function on C whose limit at z = 00 equals / . 
Therefore, 

Y{z)Y-^{z) = I or Y{z) = Y{z). 

Consider also a system {P„(z)} of orthogonal polynomials on the cross L 
with the weight ^se--^^(^), 

(E.4) / P„(z)P^(z)e-^^(^)dz = K6nm, Pn{z) = + . . . , 

where 

r ^ J y+oo ^ ,+too 
Jh ~ 27r y_oo 27r J_ 

too 

We shall assume that the nondegeneracy condition, 

•^L J jh=0,... 



(E.5) det 



7^0, for all n = 0,1,2,..., 
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is satisfied. Using tlie symmetries, z ^ —z and 2; — > one can easily show 
that (E.5) is true if, for example. 



(E.6) 



si G iR \ {0} and S3 G 



Condition (E.5) guarantees that the orthogonal set {Pn(-z)} exists. Each 
of the polynomials Pn{z) is uniquely determined from the system of linear 
equations, 



(E.7) 



I z^Pn{z)e-^^^^Uz = 0, j = 0, 1, . . . n - 1. 
7l 



Moreover, it also follows from (E.5) that 
(E.8) 

hn= [ P^{z)e-^^^^Uz = I z''Pniz)e-^^^'Uz 7^ 0, Vn = 0, 1, 2, . . . . 

Indeed, if for some n condition (E.8) is violated, the coefficients of Pn{z) would 
provide a nontrivial element of the kernel of the matrix. 



j,k=0,...,n 



Remark E.l. A truncated condition. 



(E.9) 
det 



j,k=0,...,n 



7^0, Vn = 0, 1,2, ...,m-l, m G N, 



yields the existence of the truncated orthogonal system {-Pn(-2)}™=0' 
inequalities, 

hn ^ 0, n = 0, . . . , m — 1. 

The orthogonal polynomials considered in the main text correspond to the 
choice, 



(E.IO) 



S3 = 0, si = — 27ri. 



Similar to special case (E.IO), orthogonal polynomials (E.7) satisfy the 
linear equations (cf. (1.4)-(1.6)), 



(E.ll) 
and 



zPn{z) = Pn+l{z) + RnPn-l{z), 



(E.12) Pl^{z) = NRn[t + g{Rn-l +Rn + Rn+l)]Pn-l{z) 

+ g{NRr,-2Rn-lRn)Pn-3{z), {') = 



A. 

dz 
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where 

(E.13) Rn = -^, n=l,2,.... 

The coefficients Rn satisfy the same Freud equation (1.7); i.e., 
(E.14) n = NRn[t + g{Rn-i+Rn + Rn+i)\, n = l,2,... 

with the initial conditions, 

The relation between the Riemann-Hilbert problem (E.l), (E.2) and the 
orthogonal system {Pn{z)} can be formulated as the following proposition: 

Proposition E.l (cf. [FIK2,4]). The Riemann-Hilbert problem (E.l), 
(E.2) is (uniquely) solvable for all n £ N if and only if condition (E.5) is sat- 
isfied, and hence the orthogonal system {Pn{z)} exists. The (unique) solution 
Y{z) of problem (E.l), (E.2) is given by the exact formula, 



(E.15) Y{z) = 



1 P„_i(^) -If f-.M'-""'" ,;^;' 

hn-1 " -^V / h„-i JL n-z ^/ 

where the symbol is as introduced in (E.4). 

Proof. Suppose that (E.5) takes place, and define the function Y{z) by 
the right side of equation (E.15). Then Y{z) solves (E.l), (E.2). Indeed, the 
jump equation (E.l) follows directly from the Plemelj formula for the boundary 
values of Cauchy integrals, while the orthogonality conditions, 

/ z^P„(z)e-^^Wdz = hnSnj, j = 0, 1, . . . , n, 

imply that 

^^^^-1, 0{z-') z- + 0{z—'))^ 

and hence asymptotic equation (E.2) is valid as well. 

Conversely, assume that the Riemann-Hilbert problem has a solution Y(z) 
for all n G N. The triangularity of the jump matrix S allows us to represent 
this solution in closed form (cf. [FIK4, §3.4]). In fact, asymptotic condition 
(E.2) implies that 

rFlfi^ Vt ^ fz^ + oiz^ 0{Z-^) \ 
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At the same time, the 11 and 21 components of (E.l) yield {Y-^-(z))u = 
{Y-{z))ii and {Y+{z))2i = {Y-{z))2i respectively. Using these equations and 
(E.16) we find 



(E.17) 



(y(z))ii = Pniz), {Y{z))21 = Qn-liz), 



where Pniz) and Qn-i are the polynomials of the indicated degrees with the 
only restriction that 

P„(z) = z" + ... . 

The 12 and 22 entries of equation (E.l) can be written down as the jump 
conditions, 

(y+(z))i2 - {Y-{z))u = .e-^^(^)(y_(z))n, 
(y+(z))22 - iY-{z)h2 = se-^^W(F_(z))2i, 

which together with (E.17) provide us with the following representation for the 
function Y{z): 

Qn-iMe-^^W 



(E.18) 



Y{z) 



-djjL 



\Qn-liz) -/l 

Referring again to (E.16), we conclude that the following equations, 
/ /x^P„(/x)e-^^(^)d/x = 0, j = 0, 1, . . . n - 1, 



(E.19) 
and 

(E.20) / M^Qn_i(M)e-^^(^)dM = kn-i, 



j = 0, l,...n- 1, 



hold for all n G N. 

Our next step is to obtain the nondegeneracy conditions (E.5). Suppose 

that 
(E.21) 



det 



0, for some m G {0, 1, 2, . . . , }. 



j,k=0,...,m 

Then there is a nonzero tuple, {cq, . . . , c^} such that 



Cfc / z'^z^e-^^^'Uz = 0, for all j = 0, . . . , 

k=0 



m. 



This in turn implies that the polynomials. 



Pm+liz) = Pm+liz) + CkZ^, 



k=0 
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and 

m 

Qm{z) = Qm{z) + ^ CfcZ*', 
fc=0 

satisfy (E.19), and (E.20) respectively. The pairs {Pm+i{z),Qm{z)), {Pm+i{z), 
Qm{z)), {Pm+i{z),Qm{z)), and {Pm+i{z),Qm{z)) generate, via (E.18), the 
different functions which solve the same Riemann-Hilbert problem (E.l), (E.2) 
for n = m + 1. This contradiction shows that assumption (E.21) is false. 

Condition (E.5), which has just been established, means that equations 
(E.19) and (E.20) determine polynomials Pn{z) and Q„_i(^;) uniquely for all 
n G N. Simultaneously, the equation, 

Qn-l{z) = -j^Pn-liz), 

holds. This equation and equations (E.18), (E.19) complete the proof of the 
theorem. 

Remark E.2. Let m G N. Then the solvability of the Riemann-Hilbert 
problem (E.l), (E.2) and representation (E.15) for n G {l,...,m} only need 
truncated solvability condition (E.9). 

Remark^.3. If n = 0, the Riemann-Hilbert problem (E.l), (E.2) is always 
solvable, and the solution is given by the equation, 

(E.22) ^(^)=(o r 

Prom exact representation (E.15) it follows that the solution Y{z) of the 
Riemann-Hilbert problem (E.l), (E.2) has, in fact, a complete asymptotic series 
at z = oo. 



(E.23) i + V^U"-3^ ^^oo. 




Observe also that if Y{z) is a solution of problem (E.l), (E.2), then so is the 
function, 

(-l)^asY{-z)as. 

Because of the uniqueness of the solution of the Riemann-Hilbert problem, the 
symmetry equation 

(E.24) Y{z) = {-ira3Y{-z)a3 

holds. Equation (E.24) in turn implies that in asymptotic series (E.23) all the 

even coefficients, @2i, are diagonal while all the odd coefficients, @2i+i, are off- 
diagonal. Moreover, again from (E.15), we conclude that 
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(E.25) (ei)i2 = (61)21 = 7 — , 

and hence 

(E.26) (01)12(61)21 =i?„, 

where the /i„ and i?„ are defined as in (E.8) and (E.13) respectively. 
Define 

, NV(z) 

(E.27) ^n{z) = e-^"''^TQY{z)e 



where 



To = ( Q ^-1/2 ) , and e^^" = K- 



The function ^n{z) is analytic on C \ L. It has the limits ^'^±(2;) from the it 
sides of L which satisfy the jump equation (cf. (E.l)), 



(E.28) ^n+{z) = ^n-{z)S, S = 



1 s 
1 



In virtue of (E.23), (E.25), and (E.26), the function ^n{z) has the asymptotic 
expansion. 



(E.29) ^n{z) - ( E ^ ) e-(^-"^'^^+^")'^^ z ^ 00, 

where 



\k=0 



Also, equation (E.15) yields the following representation for the function "^n^z), 

1pn{z) ^n{z) 
i'n-liz) 'fn-l{z) 



(E.31) ^n{z) 
where 



1 -NV{z) 

i'niz) = —=Pn{z)e 2 . 



and 

_ NV{p) 

Equation (E.31) provides an orthogonal polynomial representation for the 
solution of the Riemann-Hilbert problem (E.28)-(E.30). In the particular case 
(E.IO), this problem is exactly our basic Riemann-Hilbert problem (5.16)- 
(5.18), and formula (E.31) coincides with formula (5.15). Equation (5.15) 
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for the solution of the Riemann-Hilbert problem (5.16)-(5.18) was derived in 
Section 5 via the analysis of the monodromy properties of the solutions of the 
Lax pair (3.6). Matrix equations (3.6) are equivalent to the pair of the scalar 
equations (E.ll), (E.12). In the approach presented in this appendix we do not 
use these equations at all. It is the Riemann-Hilbert problem (E.l), (E.2), or 
equivalently (E.28)-(E.30), that we make now a starting point of the analysis 
of the orthogonal polynomial system (E.7). 

The Lax pair, and hence the Freud equation, for orthogonal system (E.7) 
can be obtained directly from the Riemann-Hilbert problem (E.28)-(E.30) by 
use of the general inverse monodromy problem technique of [JMU] (see also 
[FIK4]). To this end, put 

(E.32) Z{z) = e^"'^3r-i^^(^)eWWa3 = y(^)^-n^3^ 

where we denote for the sake of brevity. 

Then 

(E.33) *„(z) = e-^"'"3^oZ(z)e-'^(^)'"^ 

and by (E.23), 

(E.34) z{z)^I + ^^, z^oo. 

fe=i ^ 

The next lemma shows that the Riemann-Hilbert problem implies a pair of 
polynomial matrix differential and difference equations on ^'^(z). We shall use 
the following usual notation: If 

m 

^i^) ~ Yl ^^^^^ z^oo, 

k=—oo 

we denote by 



k=0 

the polynomial part of B(z) at infinity. 

Lemma E.l. Assume that ^niz) is a solution of the Riemann-Hilbert 
problem (E.28)-(E.30). Then ^n(-z) satisfies the polynomial 2x2 matrix differ- 
ential-difference Lax pair {cf. (3.6)), 

^ ' ' \ K{z) = NAn{z)^n{z), ^' dz' 
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with 

(E.36) An{z) = -{l/2)e-^-<T^-^o^^z{z)V\z)a^Z-\z)^ To^e^-^S 



and 
(E.37) 



1 
' ' 



where Z{z) = Zn{z) is defined as in (E.32). 

Proof. Observe that det*„(z) is an entire function, since 
det (z) = det (z) det S = det (2;) . 

In addition, 

lim det*„(z) = R-^/^ lim det Z{z) = R'^/^detl = R'^/'^. 

z — ^00 z — >co 

Hence 

det ^{z) = R-^/^ ^ 0. 
Let us check that ^{z) satisfies a matrix differential equation. Define 

Q{z) = K{z)^-\z). 
Then by (E.28) and the fact that the jump matrix S does not depend on z, 

Q+{z) = = K-{z)SS-'^-'_{z) = Q^iz), 

so that Q{z) is an entire matrix-valued function. By (E.33), 

Q{z) = e-^""3ro \z'{z)e-^^'^''^ - Z{z)W' {zjase-^^'^^' 

X e^^'^''^Z-\z)TQ^e^"''^ 
= e-^^^^To [Z'{z) - Z{z)W'{z)a^] Z''^ {z)V^^ e^'^"'' ■ 
hence Q{z) grows polynomially at infinity, and hence Q{z) is a polynomial, 

Q{z) = e-^"<^3ro| [z'{z) - Z{z)W'{z)a^] Z-\z)^ Tq ^e^"^^ 

= -{N/2)e-^»''^VQ^yZ{z)V'{z)a^Z-^{z)^ V^^e^^"^ . 
Thus we get a polynomial differential equation on ^^{z), 

= Q{z)^n{z), 
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with 



Q{z) = -{N/2)e-^"''^ro!^Y{z)V'{z)a3Y-\z)^ T^^e^-''^ 



This proves equation (E.36). Equation (E.37) is proved similarly when we 
consider the "discrete" logarithmic derivative, 



Q{z)=^n+1{Z)K'{Z), 

and take this time into account the fact that the jump matrix S does not 
depend on n. The proof is complete. 



A straightforward calculation based on asymptotic series (E.34) and sym- 
metry equation (cf. (E.24)), 

Z{z) = a3Z{-z)a3, 

leads to representations (3.7) and (3.8) for the matrices Un{z) and An{z) re- 
spectively. Also, substituting (E.33), (E.34) into the second equation in (E.35) 
and equating the terms of order z^^ , one obtains, for Freud equation 
(E.14). The latter simultaneously is a compatibility condition of the Lax pair 
(E.35) (see §5). 

One can repeat the arguments used in Section 5 for the particular case 
(E.IO) and show that the Riemann-Hilbert problem (E.29)-(E.31) represents 
the Stokes phenomenon for the matrix differential equation 

(E.38) <(z) = NAn{z)^n{z), 

^ . , / + gzRn) Rl/^[t + gz^ + g{Rn + Rn+i)] \ 

^^^> { _R]l^[t + gz^ + g{Rn-i + Rn)] f + ^ + gzRn J ' 

with Rn defined as the solution of initial problem (E.14) for the Freud equation. 

For the analysis of the general monodromy problem for equation (E.38) 
corresponding to an arbitrary solution of the Freud equation and for its relation 
to the fourth Painleve equation we refer the reader to Section 3 of [FIK2] . The 
monodromy problem corresponding to the polynomial V{z) of an arbitrary 
even degree is considered in [FIK4]. 

Remark E.4. The technique used in this appendix was first suggested in 
[FMA] for analyzing the explicit solutions of the Painleve equations. In [FIK4] 
it was applied to the case of the arbitrary even polynomial V{z). In fact, 
using the same idea one can easily reduce the analysis of an arbitrary system 
of the orthogonal polynomials {Pn{z)} on some contour L with some weight 
uj{z) to the analysis of relevant 2x2 matrix Riemann-Hilbert problem. The 
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RH problem is formulated for a 2 x 2 matrix function Y(z) which is analytic 
outside the contour L, normalized by the asymptotic condition 

1 



Y{z)z-'"'^ /, Z ^ oo, (T3 

and whose boundary values Y±{z) satisfy equation: 



-1 



This Riemann-Hilbert problem can also be used to explain the appearance 
(see, e.g., [ASM]) of the KP-typc hierarchies in the matrix models. Indeed, let 
us assume that the weight function a; (A) is of the form 

and put 

Then, the same arguments as the ones we used for proving Lemma E.l yield 
(cf. (1.12)-(1.16) in [FIK4]) the system of linear differential and difference 
equations for function ^'(2;) = 'i>{z;n,ti,t2,ts, ...), 

^{z; n + 1) = U{z)^{z; n) 
d^^[z) = A{z)^{z) 

dt,^{z)=Vk{z)^{z), fc = l,2,3,..., 

where U{z), A{z), and Vfc(z) are polynomial on z (for their exact expressions in 
terms of the corresponding i?„ see [FIK4]). The first two equations constitute 
the Lax pair for the relevant Freud equation. The compatibility conditions 
of the third equation with the different k generate the KP-type hierarchy of 
the integrable PDEs; the compatibility condition of the second and the third 
equations produces the Virasoro-type constraints; the compatibility condition 
of the first and the third equations is related to the Toda-type hierarchy and 
vertex operators. 
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